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PREFACE 


This report is divided into two parts. Part I, "Lyapunov Stability 
Theory and the Stability of Solutions to Partial Differential Equations" is 
a slight revision of the Doctoral Dissertation of Dr. Gabe R. Buis which was 
submitted to and approved by the Graduate Faculty of the School of Engineering, 
University of Pittsburgh, September, 1967. This research was directed by 
his Major Advisor, Dr. William G. Vogt, Associate Professor of Electrical 
Engineering and Principal Investigator of the grant. Part II, "Contraction 
Groups and Equivalent Norms," by William G. Vogt, Martin M. Eisen and Gabe R. 
Buis presents further extensions of some of the research reported in Part I. 

This research was supported under the National Aeronautics and 
Space Administration Grant No. NGR 39-011-039 with the University of 
Pittsburgh. 
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PART I 


LYAPUNOV STABILITY THEORY AND THE STABILITY OF SOLUTIONS 
TO PARTIAL DIFFERENTIAL EQUATIONS 


by 


Gabe Rinse Buis 

• Department of Electrical Engineering 
University of Pittsburgh 


ABSTRACT 


Lyapunov stability theory is generalized to semi-groups and groups of 
linear operators in Hilbert spaces. The existence of a Lyapunov functional 
is sufficient for the asymptotic stability of semi-groups and necessary and 
sufficient for the asymptotic stability of groups. This theory is applied 
to a class of partial differential equations, yielding stability conditions 
which are mathematically rigorous. 
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SUMMARY 


Lyapunov’s stability theory has become very important in the stability 
analysis of solutions to ordinary differential equations. Its extension to 
partial differential equations has been restricted to a few scattered 
applications which generally lack mathematical rigor. The complex mathemat- 
ical nature of partial differential equations makes extension of the stability 
theory to partial differential equations very difficult. The objective of 
this dissertation is to develop a mathematically rigorous Lyapunov stability 
theory for a class of partial differential equations. 

The approach taken to this problem is that of generalizing Lyapunov's 
Direct Method for ordinary differential equations to a class of operator 
differential equations. This generalization is based on the fundamental 
solution structure of groups and semi-groups, which is not restricted to 
ordinary differential equations. 

This leads to the formulation of sufficient conditions on an operator 
to generate stable or asymptotically stable semi-groups and necessary and 
sufficient conditions on an operator to generate asymptotically stable groups. 

The semi-group and group structures enable one to associate with the 
infinitesimal generator an operator differential equation. The derived 
stability conditions extend then to the stability of the null solution of 
this differential equation. The generators of these semi-groups and groups 
constitute a large class of bounded and unbounded operators. 

The second stage is to formulate partial differential equations in 
the framework of the operator differential equations. Starting with a formal 
partial differential operator, distributions are introduced to define the 
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extended operator in a complete space. The next step is to associate boundary 
conditions with the formal partial differential operator. This is limited to 
the so-called Dirichlet boundary conditions, which are important for many 
stability investigations. 

For the class of strongly elliptic partial differential operators the 
domain and range are formulated in terms of Sobolev spaces. Subsequently, the 
developed stability theory is applied to a class of evolution equations with 
strongly elliptic partial differential operators and Dirichlet boundary 
conditions, giving sufficient conditions for asymptotic stability of the null 
solution. 

A similar formulation is given for a class of wave equations. However, 

for this case the necessary and sufficient conditions for asymptotic stability 

of the null solution are established by developing a suitable Lyapunov 

Functional. At the same time this proves that the solutions possess the group 

* 

property. Various applications are given. The stability analysis is further 
extended to some nonlinear partial differential equations. With the exception 
of these last results the emphasis has been on a mathematically rigorous 
formulation of the stability problem. 



I. INTRODUCTION 


In recent years Lyapunov stability theory has become an important tool 

in the stability analysis of solutions to linear and nonlinear ordinary 

( 1 )* 

differential equations. The original work of Lyapunov has generated many 
contributions to the stability theory of solutions to ordinary differential 

equations and has provided many applications. Two of the more significant ones 

( 2 ) ( 3 ) 

are the paper by Kalman and Bertram and the book by LaSalle and Lefschetz' ; 

Where these contributions involve Lyapunov's Direct Method, the central problem 

becomes the construction of a Lyapunov function. For nonlinear systems this is 

generally very difficult. 

Although the development of Lyapunov stability theory and applications 
to solutions of ordinary differential equations has progressed rapidly, its 
application to solutions of partial differential equations has remained 
limited. This is the case despite the fact that many physical systems must be 
represented by partial differential equations. However, the difficulties 
encountered in applying Lyapunov stability theory to partial differential 
equations parallel those in establishing the existence and uniqueness 
properties of solutions to partial differential equations. 

Yet stability remains one of the most important properties of distributed 
parameter systems. Many of the stability results for partial differential 
equations are obtained by using methods of approximation. These methods might 
not give sufficient conditions for stability except in the case of infinites- 
imally small perturbations. Stability can be defined in many different ways. 

In reference to stability theory it will be interpreted here as stability in 

*Parenthetical references placed superior to the line of text 
refer to the bibliography. 
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the sense of Lyapunov: a system is said to be stable if for sufficiently 

small initial perturbations the solutions remain close to the original solution 
for all future time. The advantages of Lyapunov's Direct Method over 
approximate methods are many. Most important, Lyapunov's Direct Method employs 
the system equations directly without resorting to approximations. The method 
also allows mathematical rigor and in principle the introduction of nonlinear- 
ities. 

On the other hand, there is the importance of partial differential 
equations in the fields of reactor physics, hydrodynamics, magnetohydrodynaraics , 
control processes, etc. These facts certainly motivate an investigation of 
possible ways to extend Lyapunov stability theory for the stability analysis 
of solutions to partial differential equations. The following section gives 
a survey of the significant results obtained so far. This survey shows how 
limited these results are. The survey also points out the lack of mathematical 
rigor in many of the applications. This section is followed by a short out- 
line of general problem areas in the stability study of solutions to partial 
differential equations. 

A. Review of the Literature 

Many stability results for partial differential equations are derived 
by using approximate methods. The basis for these approximate methods is the 
reduction of the partial differential equations to a system of ordinary 
differential equations. This can be done by either approximating the model by 
one having a finite number of degrees of freedom via spatial discretization or 
by assuming a harmonic time dependence. The first case allows the application 
of the well-known techniques for analyzing the stability of ordinary differen- 
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tial equations; in particular, for infinitesimally small perturbations, which 
is presented as a justification for the system linearization. 

In the second case a modal analysis is in general necessary. To 
achieve this, use is made of the Galerkin process which is based on a trunca- 
tion of the modal expansion. A linearization again limits the amount of work 
involved. The use of these methods is wide-spread and well published. Since 
this approach does not constitute the subject of this thesis, reference will 
just be made to the works by Bolotin and Eckhaus^ as the most recently 

published books. 

As distinct from these approximate methods, Lyapunov's Direct Method 
deals directly with the system of partial differential equations without 
resorting to approximation. Moreover, it is potentially applicable for the 
stability analysis of nonlinear systems. Thus it is not surprising that 
attmepts have been made to apply Lyapunov's Direct Method to derive sufficient 
conditions for the stability of equilibrium solutions of systems of partial 
differential equations. A step toward applying Lyapunov's Direct Method to 
partial differential equations was made by Massera^ who extended this method 
to denumerably infinite systmes of ordinary differential equations. 

The application of Lyapunov's Direct Method for the stability analysis 
of solutions to partial differential equations requires a generalization of 
the method to function spaces in which a metric p is defined. Consequently, 
the concepts of stability must be defined in terms of this metric. A general 
stability theory now based on the existence of a Lyapunov functional is 

ro\ 

established by Zubov v J for the invariant sets of dynamical systems iii general 
metric spaces. Zubov employs this theory in Chapter 5 to derive results for 
systems of partial differential equations. 


6 



The most general type of system that has been considered is of the form 


3 jj(t,x) 

= Lu(t,x) (1-1) 

where iu(t,x) is an n-dimensional vector valued function defined over some 
region ft of an m-dimensional Euclidean space E m . L is a linear or nonlinear 
matrix differential operation defined on ft. To specify solutions, a set of 
boundary conditions must be given. In addition, a solution will depend on 
some initial function Uq(x) belonging to an n-dimensional space of initial 
functions . 

The application of Lyapunov stability theory for the determination of 
stability conditions for equilibrium solutions of (1-1) is almost entirely 

(Q\ 

based on the work of Zubov v . However, in all those cases the validity of 
the results depends on the system being a dynamical system, i.e., on the fact 
that the solutions possess the group property, or to a somewhat lesser extent, 
the semi-group property (t > 0 only). The mathematical justification of this 

fact is either extremely vague or omitted. 

( 8 } 

Zubov'' ' establishes quite conclusive results for systems of the form 
3u(t,x) 3u 

T t ' ’ - Si 5 (I ‘ 2) 

the solutions of which, constitute, under suitable assumptions, a dynamical 
system. Brayton and Miranker v 1 apply his results to establish stability 
conditions for a nonlinear system representing an electrical circuit, without 
properly verifying the exact conditions for a dynamical system. Blodgett 
takes a chemical reactor model to apply Zubov's results. 

Zubov also compares the stability properties of the trivial solutions 
of the system 
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(1-3) 


3 u 

ft 


s 


£ e (u 1> 


k 3 u 

..., u n ) + I b i (s=l,2,...,n) 
i=l i 


and the related system of ordinary differential equations 
du 

■jjr- = f g u n ) (s=l,2,...,n) (1-4) 

He shows that the asymptotic stability of the trivial solution of (1-4) assures 
the asymptotic stability of the trivial solution of (1-3). A similar result 
relates the stability behavior of the equilibrium of the system of partial 
differential equations of higher order 


3 u M 
3 1 m 

I v° 

j=l 3 


+. . .+ a 


in u 


S. *” ... 9 : m 

X X 

m 


d-5) 


to the stability of the equilibrium of the system 

d ti 

~3t ’ Ai- 


d-6) 


The nature of the derived results is very theoretical and often 

( 11 ) 

difficult to implement in practical applications. Hsu v ' applied Zubov's 

results to a nuclear reactor system, but does not verify the dynamical system 

( 12 ) 

properties. Wang in a kind of survey paper, considers those operators 

which are infinitesimal generators of semi-groups. However, in this case, 

the conditions are only sufficient for stability. 

Another class of systems frequently encountered is of the form 
2 

3 u(t,x) 3u(t,x) 

5 + v + L u(t,x) = £ (1-7) 


3 t‘ 


3 t 


with ij(t,3t) and ]L as under (1-1) . Although (1-7) can be reduced to the form 
(1-1), it has the distinct advantage that the Lyapunov functional can readily 
be derived from the total system energy, again giving only sufficient condi- 
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tions for stability. Most other contributions are not as general as the ones 
above, but reflect more direct applications to specific problems. 


Movchan 


(13) 


considered the equation 

2 4 2 

3jui + 3_u _ a 9_u = Q 

9 t 2 9 x 4 9 x 2 


(1-8) 


with the boundary conditions: 


3 2 

u = — r- = 0 for x=0 and x=l. 

9 x^ 


By defining the metric p in a suitable manner he is able to verify results from 
the theory of vibrations of plates by taking as Lyapunov functional: 


V(u) = / (u + a u + u ) d x. 

^ xx x t 


,(14) 


Similarly Movchan v verifies classical stability results for a system of 
hinged rectangular plates under compression, the deflection of which, u(x,y,t), 
is given by the dimensionless equation: 

4 


2 4 

LJL + a LJL+ 2 IH 

21 A z 2 2 

9 t 9 x^ 9 x 9 y 

with boundary conditions: 


4 2 2 

, 9 u 2 9 u 2 9 u n 

+ a 2 ~ 4-’ a 3-2' “4 7^ = 

9 y 9 x 9 y 


(1-9) 


(15) 


2 2 

u = — -r = 0 at x=0, x=l, u = — ^ = 0 at y = 0, y = 1. 
9 x X 9 y^ 


Wang SJ "' / studies the stability of a simplified flexible vehicle with a 
dimensionless equation of perturbed motion about its equilibrium state as 
given by 


m 


(x)v 2 £ 2 ~ u % - x) + v A.U.x) - u - ( k ^l 


9 t 


9 t 


(1-10) 


+ ~ 2 (ei(x) ~ u( - 2 - > - x) ) = 0 

9 x 9 x 


and boundary conditions 


u(t,0) = 0, I = 0 ; EI(x) 8 — U< - t? - X - I = 0 

9 x *x=0 9 -k ' x=l 
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But he does not consider the existence problems of the solutions involved. 

Parks' 1 1 applies Lyapunov's Direct Method to the panel flutter 
problem. The equation in dimensionless form is given by 

„ + is + d iia - f l-s + h — = o (i-u) 

at at ax ax ax 

and boundary conditions 

~2 

u = — = 0 for x = 0 and x =1. 

3 x 

In these applications an important role in deriving the sufficient 

conditions for stability is played by integral inequalities. The results 

obtained are presented without mathematical rigor. Another deficiency is the 

general interpretation of the conditions obtained for stability. 

In order to discuss stability in a meaningful sense it is often 

necessary to put restrictions on the initial states. Although Volkov^^ 

implied this in an earlier work, the idea of introducing a second metric for 

this purpose seems to have been originated by Movchan' . Stability is then 

(19) 

defined in terms of the two metrics, rather than one. Wang uses this 

concept in a stability analysis of elastic and aeroelastic systems. 

The work of Lakshmikanthara^^ is closely related to Zubov's results. 

( 21 ) 

A practical application can be found in Wei's paper' ' in which the stability 
of a system of partial differential equations describing the first order 
chemical reaction in the presence of a catalyst Is analyzed. This system can 


be reduced to a pair of identical partial differential equations of the form 


9 u 
9 t 


9 2 u 2 
— r- - <J> u exp 

9x 


gu(l-u) 


1 + B(l-u) 


(1-12) 
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with 0 < x < 1 and boundary conditions: 

|* I = 0 and u(t,l) = 1. 

'x=0 

After linearizing (1-12) the Euclidean metric has been taken as a Lyapunov 
functional. 

This survey of all the significant contributions indicates that many 
problems concerning the application of Lyapunov stability theory to partial 
differential equations remain unsolved. In the next section the problems in- 
volved will be analyzed in more detail, and the most logical approach to 
further applications will be established. 

B. General Problem Areas 

The literature survey shows that the main theoretical contribution to 
Lyapunov stability theory for solutions to partial differential equations is 
given by Zubov' 1 1 . The important restriction on his result is the requirement 
that the system of partial differential equations defines a dynamical system. 
This implies that L of (1-1) must generate an operator which possesses the 
group property. However, many partial differential operators can 
only be defined as the infinitesimal generators of semi-groups. Thus the 
"group property" will only be satisfied for t > 0. Since stability is usually 
concerned with the properties of positive half trajectories, it seems natural 
that a Lyapunov stability theory can be formulated for systems having only the 
semi-group property. This suggests the development of a general stability 
theory for the class of operators generating a semi-group. This must be 
followed by the formulation of L in terms of this class of operators and where 
possible the results so obtained should be related to the actual boundary 
value problem. 
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After solving the problem for the time-invariant case there arises the 
possibility of extending the results to the time- invariant case. 

Another extension can be directed towards the development of a 
Lyapunov stability theory for nonlinear partial differential equations. The 
formulation of a Lyapunov stability theory for semi-groups appears a natural 
start for research in this field. This is because the infinitesimal generators 
of contraction semi-groups, which are stable, are the linear dissipative opera- 
tors. The natural extension of these linear dissipative operators for the 

( 22 ) 

nonlinear case are the monotone operators • The question arises, therefore, 
if it is possible to extend the Lyapunov stability theory for systems with 
dissipative operators to include monotone operators as well. 

Many of the above problems touch on research areas in mathematics which 
are currently being explored. It is expected that many of the above problems 
will provide an emerging field of future research in stability theory. 
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II. STATEMENT OF THE PROBLEM 


A. Formal Partial Differential Equations 

Many of the physical problems that are formulated by partial differen- 
tial equations can be described formally by the general equation 

3 u (t , x) 

-ft h. x £ ft (II-l) 

where u(t,.x) is an n-vector function and L is a matrix whose elements are 
linear or nonlinear differential operators specified on a bounded connected 
open subset ft of an m-dimensional Euclidean space, E m . The parameters of 
can be space and time dependent. 

In order to specify solutions to (II— 1) a set of additional constraints 
or boundary conditions must be given by 

Hu(t,x’) = 0 x' e 3ft (II-2) 

where H is a matrix whose elements are specified differential operators and 
3 ft is the boundary of ft, = ft + 3 ft. 

In addition to the boundary conditions, solutions to (II-l) will 
depend on some initial functions ( x) . It will be assumed for the moment 

that, given some initial function (x.) belonging to some normed linear space 
H, it can be shown that solutions to (II-l) and (II-2) exist and belong to H. 

A solution to (II-l) and (II-2) will be designated as u(t,x;uQOc) ,tg) , that is, 
the solution starting at t Q and with initial conditions Uq(x), u(tg »x>_Uq (x) , 
Uq(x). 

The solutions that are of particular interest in stability studies are 
the equilibrium solutions, u (x). The equilibrium solutions can be defined 

"'6CJ 

as: 
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Definition II-l . An equilibrium solution, u (x) is a solution of (II-l) and 
9u (t,x) 

(II-2) such that — \ = 0 for all t > t n and all xeft. Thus 

, to) = u c ^ (x) for all t > t^. This is the same as determining the 

u (x) such that L u (x) = 0 and H u (x f ) = 0, x f e3ft. 

Stability can be defined in many different ways; however, stability 
will be defined here in the sense of Lyapunov. 

Definition II-2. The equilibrium solution u_(x) of (II-l) and (II-2) is said 
" ec l 

to be stable in the sense of Lyapunov if for every real number e > 0, there 

exists a real number 6 > 0 such that for u 0 (x)eH ,] lu^Cx)-^ (x) I I < 6 implies 

— "U — — u — — eq — 

| |u(t,x;uQ(x),t 0 ) - UgqCx)! | < e for all t > t Q . 

Definition II-3 . The equilibrium solution u^Cx) of (II-l) and (II-2) is said 

to be asymptotically stable if it is stable and in addition | lj£(t>2£>UQ(x) £ q) - 

u (x)||-*0 as t-x». 

-eq 


It follows from the literature survey that many authors have investiga- 
ted problems similar to the one formulated above using Lyapunov’s Direct Method. 
In many cases no attention has been paid to the question of existence of solu- 
tions to (II-l) and (II-2). Furthermore Lyapunov’s Direct Method for ordinary 
differential equations is based on the properties of bounded operators in 
finite dimensional space. In solving the stability problem for (II-l) - (II-2), 
it is generally assumed that the functions are defined on the proper spaces, 
i.e. , the functions possess enough continuous derivatives. It is furthermore 
assumed that all operations of differentiation, integration, etc., can formally 
be carried out and that the functions resulting from these operations still 
belong to the defined normed linear space, i.e., that the space is complete. 

In order to satisfy all these requirements and thus conduct a rigorous 
investigation of the existence and stability properties of solutions to (II-l) 
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and (II-2), the formulation of the problem must be carefully scrutinized. This 
investigation makes extensive use of the concepts of functional analysis, in 
particular, the analysis of functions defined in complete spaces. For this 
purpose the formal partial differential equation (II-l), that is, without 
boundary conditions, is considered and the properties of the formal partial 
differential operator are studied in spaces that possess the necessary 
differentiability properties. As a part of the required technique, it is 
necessary to complete the spaces of , differentiable functions. For this purpose 
it is necessary to introduce suitably defined generalized derivatives. 

Once this process is carried out (II-l) can be interpreted as a general 
operator differential equation. A Lyapunov stability theory will be developed 
for a certain class of operator differential equations, which is analogous to 
Lyapunov's Direct Method for ordinary differential equations. The formal 
partial differential operator is then a restriction of a differential operator 
from this class. Finally there remains the problem of relating the stability 
properties of the formal partial differential equation to those of the original 
boundary value problem. In the next section the general operator differential 
equation is evaluated. 

B. Evaluation to Operator Differential Equations 

Lyapunov's Direct Method for finite dimensional systems of ordinary 
differential equations is based on certain fundamental properties of the 
solutions, which all.systems of finite dimensional ordinary differential 
equations have in common. Similar properties hold for certain classes of 
operator differential equations to which certain classes of formal partial 
differential equations belong. 
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For the characterization of the operator differential equation 

consider 

d_x 

- A x (II-3) 

Let (II-3) be valid for all jceX where X is some n-dimensional Euclidean space, 
E n , and A a linear operator on X into X. Thus, A is a bounded linear operator. 
Then (II-3) becomes a stationary linear ordinary differential equation for 
which a matrix representation is obtained by choosing a basis in E n . A 
Lyapunov stability theory has been developed for these equations. This theory 
is based on the properties of bounded operators in finite dimensional spaces. 

If X is a general normed linear space, then A no longer need have these 
properties; in fact, A might be unbounded and (II-3) should be specified as 

dx 

^=Ax (xef(A) C X; A : P (A) -*■ R(A)C X) (II-4) 

where P(A) is the domain of A and R(A) is the range of A. For these unbounded 
operators on general function spaces, a mathematical theory has emerged in which 
the properties of the solutions are studied on one-parameter families of 
bounded linear operators, the semi-groups and groups. 

Lyapunov's Direct Method for ordinary differential equations requires 
the construction of a Lyapunov function . For a general operator differential 
equation (II-4) this requires the construction of a Lyapunov functional . The 
stability properties follow from an evaluation of the time derivative of the 
Lyapunov functional along the solutions. For the Lyapunov functional the 
time derivative must formally be defined, and there is a possibility that it 
does not exist for all jceX. In the context of semi-groups this time derivative 
can easily be derived. 
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In the following it should be clear from the context whether x must be 
interpreted as an element of the normed linear space X (in the general 
theoretical development) or whether x is a space variable in E m , the space on 
which the differential operator is defined. 

A large class of formal linear stationary partial differential equations 

3u(t,x) 

9“t = —(2.) Ji (tjX.) x e ft (II-5) 

where ft is a bounded connected open subset of an m-dimensional Euclidean space, 
E m , and Mjx) is a formal partial differential operator on the space variable 
x, is subsummed in a family of general operator differential equations (II-4) 
which generate a semi-group or a group. Thus, the operator A in (II-4) is an 
extension of L^jc) in (II-5) and coincides with _L(x) when the functions u are 
sufficiently smooth. 

The main problem is to deduce from (II-5) the form of the space X, 
the linear operator A, the domain V (A) and the range R (A) . This in itself is 
not always easy. 

C. A Stability Theory 

In the preceding sections the difficulties of establishing a Lyapunov 
stability theory for partial differential equations (II-l) and (II-2) have 
been discussed. On the other hand, the general linear operator differential 
equation (II-4) can be considered as a generalization of a system of linear 
ordinary differential equations for which a Lyapunov stability theory exists. 

In order to derive a Lyapunov stability theory which can rigorously be applied 
to a class of partial differential equations, the stability problem will be 
formulated as follows: 
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1. Develop a Lyapunov stability theory for a class of operator 
differential equations (II-4) similar to Lyapunov's Direct Method for ordinary 
differential equations. 

2. Extend the formal partial differential equations to operator 
equations on complete function spaces in such a way that all mathematical 
operations can be carried out rigorously. 

3. Associate with the formal partial differential operator, so defined, 
a boundary value problem and give a formulation of the problem in terms of the 
general operator differential equation for which the stability theory is 
developed. 

4. Give applications of the developed stability theory to specific 
boundary value problems. 

Even though the scope of these problems is limited to a small class of 
partial differential equations, which are linear and stationary, they 
constitute one of the first developments of a mathematically rigorous approach 
to the stability investigations of partial differential equations. 

D. Contributions to the Problem 

The research concerning the stability properties of the solutions to 
certain classes of partial differential equations was carried out as part of 
a research project sponsored by the National Aeronautics and Space Administra- 
tion under Grant Number NGR 39-011-039 with the University of Pittsburgh. 

Dr. William G. Vogt, Associate Professor of Electrical Engineering and 
Principal Investigator for this project originated this research effort and 
in the course of the research has contributed considerably to the results 
reported in this thesis. 
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The approach to the stability problem and the results obtained were 
developed in the course of numerous discussions between Dr. Vogt and the 
author. Since the results reported are to a large extent obtained as a joint 
effort, it is difficult to isolate many single results as principally Dr. Vogt’s 
or the author’s. Some of these are indicated below. 

Most of the research carried out so far on this subject lacked a solid 
mathematical justification. The initial efforts by the author gave a 
verification, using Lyapunov stability theory, of stability results obtained 
by Eckhaus ^ who used approximate methods. These results were eventually 
established rigorously for the linear case (the case for which the results of 
Eckhaus are certain to hold) and are given in Chapter IX. 

The important contribution at this stage was the introduction of the 
concept of equivalent inner products by the author. However, its use was 
primarily aimed at obtaining self-adjointness properties for the formal 
operator involved and improving the estimates obtained by the use of integral 
inequalities. The abstraction of this notion and its final implementation in 
a stability theory for dissipative operators is largely due to Dr. Vogt. 

Lyapunov’s Direct Method for systems of ordinary differential equations 
is based on some fundamental properties, i.e., group properties, of the solu- 
tions. A first observation showed that this group property is not limited to 
finite dimensional systems of ordinary differential equations, but is common 
to a larger class of systems. Moreover, a still larger class of systems 
possesses the more general semi-group property. This led to a formulation of 
a Lyapunov stability theory for groups and semi-groups. The stability 
properties are in this case directly coupled to the trajectory structure, i.e., 
the solutions. The main goal became to express the conditions for the 
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stability of the solutions in terms of the operator generating the group or 
semi-group — in particular, those generating contraction groups or semi-groups. 

A large class of such operators is formed by the bounded and unbounded 
dissipative operators. 

Dissipativity, in a particular case, is an inner product property while 
stability is a norm property. This implies a Hilbert space theory. Since the 
stability properties are invariant under equivalent norms, the introduction of 
the principle of equivalent inner products mentioned above is crucial to link 
the dissipativity property to the stability properties. Two inner products 
are equivalent if and only if their induced norms are equivalent. 

This subsequently allows the formulation of sufficient conditions on 
an operator to generate stable or asymptotically stable serai-groups and 
necessary and sufficient conditions on an operator to generate stable or 
asymptotically stable groups. Dr. Vogt proved an important norm property of 
groups which is used in the proof of the last statement. 

The semi-group (group) structure enables one to formulate an operator 
differential equation with its infinitesimal generator. The derived stability 
conditions extend to the stability of the null solutions of this differential 
equation. The stability theorems are given in Chapter V. 

The conditions to be imposed on the operator do not restrict the class 
of operators to bounded operators. Thus a stability theory is established for 
a large class of operator differential equations defined in a Hilbert space. 

The second stage is to formulate partial differential equations in the 
framework °f tbe operator differential equations. A kind of synthesis proced- 
ure is used here. 
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Starting out with a formal partial differential operator, the first 
requirement is that the extended operator be defined in a complete space. This 
is done by introducing distributions. The next step is to associate boundary 
conditions with the formal partial differential operator. In the context of 
the thesis this is limited to the so-called Dirichlet boundary conditions, 
which are important for many stability investigations. 

The class of formal partial differential operators is limited to the 
strongly elliptic partial differential operators. The domain and range of the 
extended operator are subsequently formulated in terms of Sobolev spaces. In 
Chapter VII, the stability theory is applied to a class of evolution equations 
with strongly elliptic partial differential operators and Dirichlet boundary 
conditions. 

A similar formulation of a class of wave equations is presented in 
Chapter VIII. The development of a Lyapunov functional for this class of wave 
equations is considered one of the main contributions of this thesis. It 
establishes conditions for asymptotic stability of the null solution and the 
group property of the solutions to the wave equations. 

Chapter IX is devoted to a number of applications which illustrate 
various aspects of the developed stability theory. With the exception of 
Chapter X, almost entirely the author's work, where some nonlinear systems are 
formally analyzed, the emphasis has been on a mathematically rigorous 
formulation of the stability problem. Even though only the stability proper- 
ties for a class of linear, stationary, partial differential equations are 
established here, the directions to be pursued for enlarging this class have 
been opened to further research. These results appear to be a significant 
contribution toward a rigorous Lyapunov stability theory for a more general 
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class of partial differential equations including nonlinear partial 
differential equations. 
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III. MATHEMATICAL PRELIMINARIES 


Engineers have become familiar with the concepts of linear vector 
spaces and linear transformations on these spaces through the introduction of 
the state space approach in control theory and the development of a Lyapunov 
stability theory for ordinary differential equations. The stability study of 
general operator differential equations adds another dimension to these 
concepts. Functional analysis is a fundamental tool in the study of operators 
defined in general function spaces. In the following sections some of the 
basic notions and properties that are important in the stability analysis of 
operator differential equations will be introduced. A more detailed treatment 
and examples can be found in any book on functional analysis, for example. 
Reference 23. 

A. Normed Linear Spaces 

In the previous chapter the general operator differential equation was 
introduced by defining A as an operator in a normed linear space. Such a 
space is defined as follows: 

Definition III-l . Let X be a vector space over the field of real or complex 
numbers. A norm on X, denoted by | | • | | , is a real-valued function on X with 
the following properties: 

i. | |x| | > 0 for all xeX. 

ii. X f £ implies | |x| | ^ 0. 

iii. | | ax | | = |a| | |x| |, a some real or complex scalar, 

iv. | |x + zl 1 = i 12E.I I + I Ixl I (triangle inequality). 

The vector space X, together with a norm on X, is called a normed linear space. 
When the scalars over X are the reals, X is called a real normed linear space. 
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The finite dimensional real Euclidean space, R n , can be made into a 
normed .linear space by defining the norm by 

Mill - ([ xj) 1/2 V x e R n . 

i=l 1 

The norm for a space can be defined in more than one way. Thus another norm 
for R n is defined by 

| |xj | = sup |x ± | V x e R n . 

These two norms create two different normed linear spaces, with possibly 
different basic properties. The introduction of a norm for R n makes (II-3) a 
special case of the general linear operator differential equations as defined 
by (II-4) . 

The norm of X induces a metric or distance d which is defined by 

d (x,£) = [|x - Z M x, y eX 

i.e., the distance between two elements jc and ^.eX is given by | |i “ ill * 
According to their properties one can distinguish different classes of normed 
linear spaces. First of all the pre-Hilbert spaces are defined: 

Definitions III-2 . A real or complex normed linear space X is called a 
pre-Hilbert space if its norm satisfies the condition 

lli + zll 2 + ! Is - 2.1 1 2 » 2(11*1 1 2 + Mill 2 ). 

The more important normed linear spaces are the complete normed linear 

spaces or Banach spaces. A Banach space X is a normed linear space in which 

every Cauchy sequence converges with respect to the norm to a limit point in 

X. A complete pre-Hilbert space is called a Hilbert space. The norm in the 

1 /2 

Hilbert space is the one induced by an inner product, | |xj | = <Z»x> • The 

inner product, <•»*>, is defined in a real pre-Hilbert space by: 

< 2£»z > ■ j Mlz + zM 2 - II* - zM 2 )* 
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and in a complex pre-Hilbert space by: 

<x,p = < 2 i»X > i + i< Z» i Z > i 

where < Z»Z > i = J ( I I* + zl | 2 - Mz ~ zl | 2 ) • 

This inner product has the following properties: 
i. < a jc,^> = a •cxjjp ' a real or complex 

ii. < jc + £,z> = < x> z> + < y,z> 

iii. < = < Z’Z^* f° r rea ^ space (= < y_,x> for complex space) 

iv. < jc,x> > ® whenever x ^ 0. 

The finite dimensional real Euclidean space R n , is made into a Hilbert space 
by defining the inner product by the following finite series: 

< Z»Z> = Vl + X 2 y 2 + ••• + in- 
stability is defined in terms of the norm chosen for the linear space. It 

might be expected that the stability properties will depend on the particular 

norm selected for the space. The norm can in general be chosen in many 

different ways. The stability properties are preserved from one space to 

another if the norms are equivalent. Two norms, ||*||^ and M'l^* are 

equivalent if there exist constants c^ and c^, °°>C2>c^>0 such that 

c ! I Izl I ± < Mil 1 2 = c 2 1 l-l 1 1 ( v i eX ) (HI- 1 ) 

Moreover the spaces and X2 denoted by X^ = (X; | |. | |^) and X2 = (X; | |. | 
respectively are topologically equivalent. 

For systems described on finite dimensional Euclidean spaces all 
norms are equivalent. Thus stability with respect to one norm implies stabil- 
ity with respect to another, i.e., stability in one normed space implies 
stability in another. The stability problem evolves basically around the 
selection of a norm, i.e., the normed linear space. For Hilbert spaces 
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this reduces to the selection of an inner product. 

Since the stability theory to be developed concerns primarily 

dissipative operators , which are defined in terms of the inner product, 

Hilbert spaces are very important. Dissipative operators can be studied in 

somewhat more general spaces. These so-called semi- inner product spaces were 

( 25 ^ (26} 

introduced by Lumer v . Lumer and Phillips' 1 J studied dissipative operators 
in these particular spaces. 

Definition III-3 . A semi-inner product space is defined on a complex or real 
vector space X with norm | | . | | if to each pair x,£ e X there corresponds a 
complex or real number [_x,^_] such that 

i. [x + y_,_z] = [x,z.] + [£,zj 

ii. [ax,y] = a[x,^l 

iii. [x,x] = | |x| | 2 for _x ^ _0 

iv. |tx,zll 2 = 

Any Banach space can be made into a semi-inner product space. In particular, 
for a Hilbert space the only semi- inner product is the usual inner product. 

The relevant properties of normed linear spaces can be found in any 
book on functional analysis and the papers given as references. Where 
necessary in the following chapters, these properties will be recalled. In 
the next section, linear operators and their properties in these various 
spaces are introduced. 

B. Linear Operators 

(23) (27) 

This section is devoted to the characterization of operators * . 

Let X and Y be two linear spaces over the same real or complex field. The 
objective is to characterize the mappings ^T: 2L Z where x. £ X and ^ e Y, 
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Definition III-4 « The mapping TF: x = T(x) defined on a linear manifold D 
of X and taking values in Y with the property T(ct x^ + B x 2 ) = aCTx^) + B (Txp) 
is called a linear operator on D C X into Y. 

D = P(T), the set in X on which T? operates, is called the domain of T. 
The set in Y which results from the operation of on V(T) C X is called the 
range of R(£), 

R(T) = {^eY?^=Tx, xs V(T)}. 

In certain applications the null space W(T) is used as defined by 

W(T) = {x e P(T) 5 T x = 0}. 

If the range R(jC) is contained in the scalar field K, then T^ is called a 
linear functional. If a linear operator T_ gives a one-to-one map of V(T) onto 
R(T), then the inverse map ]T -1 gives a linear operator on R(T) onto P(T): 

T^ T^ jc = x for x e V(T) and T - '*' ^ = y for e R(T) . 
is the inverse of T. 

Proposition III-l . A linear operator T_ admits the inverse T^ if and only if 
TC jc = 0 implies x - 0. 

Definition III-5 . Let T^ and be linear operators with domains P(T^) and 

V(T 2 ) both contained in a linear space X, and ranges R(T^) and RCT^) both 
contained in a linear space Y. Then T^ = T? 2 if and only if P(T^) = PCT^) and 

T^ x = T 2 x for all x e ^(T^) = PCr 2 ). If and ^l— = ^2— for a11 

x e V (T^) , then is called an extension of T^ and T^ a restriction of Tg, 

written as T^ C 

An important role is played by the bounded linear operators. Let X 
and Y be normed linear spaces . A linear operator T^ with domain in X and 
range in Y is bounded if there exists a positive constant M < » such that for 
all x e V (T) 
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'C is continuous at a point x e X and T_ is uniformly continuous in X are 
equivalent statements. The class of all bounded linear operators on X into Y 
is designated as L(X,Y) i.e., T eL(X,Y) P(T) = X, R(T) C Y, T is bounded. 
Proposition III-2 . Let X and Y be normed linear spaces. Then a linear 
operator on P(T) C X into Y admits a continuous inverse T^ if and only if 
there exists a positive constant y such that 

I I 1 2il I = Y I l x l I f° r every x e P(T) . 

Definition III-6 . If T is a bounded linear operator on a normed linear space 
X into a normed linear space Y, then its norm is defined by 

I 111 I = sup ||T x| | - sup | |T x | | 

iiiii- 1 s* Turn 

An extension of the notion of a bounded linear operator is that of a 
closed linear operator. The definition is based on the notion of graph of 
T, G(T) . 

Definition III-7 . The product space X x Y is defined as the normed linear 
space of all ordered pairs (x,jyO , x e X, ^ e Y, with the usual definitions of 
addition and scalar multiplication and with norm given by 

IIQe»z>I I = max { l l x l I* I lxl I 

Definition III-8 . The graph GOT) of T_ is the set {(x t Tx) |x e V ( T) } . Since 
1C is linear, G(1C) is a subspace of X x Y. If the graph of T^ is closed in 
X x Y, then 1C is said to be closed in X. When there is no ambiguity concern- 
ing the space X, T^ is said to be closed. 

The following remarks can be made: 

i. T^ is closed if and only if {x } in P(T), x -*■ x, 1 x -»■ jr, imply 


ii. If T is 1-1 and closed, then T_**^ is closed, 

iii. The null space of a closed operator is closed, 

iv. If V(T) is closed and £ is continuous, then T^ is closed, 

v. The continuity of T^ does not necessarily imply that T^ is closed. 

T is closed does not necessarily imply that T_ is continuous. 
Closed-Graph Theorem . A closed linear operator mapping a Banach space into a 
Banach space is continuous. 

Some additional properties of closed operators are introduced after 
defining the adjoint operator. 

For a Banach space X the conjugate or dual space, the Banach space of 

bounded linear functionals on X, is indicated by X' . Let T_ be a linear 

operator mapping the normed linear space X into the normed linear space Y and 

with domain dense in X. Then the conjugate of ]T is denoted by T' . The 

(23) 

definition can be found in . 

If X and Y are Hilbert spaces, then the notion of conjugate operator 
of T^ can be extended to that of adjoint operator of JF which is denoted by T*. 
The operators one generally deals with are mappings from subsets of a Hilbert 
space X into X.. Let T^ be such a linear operator. The adjoint T* of T with 
respect to X is defined by 

< T x, j> = < x,T* Z > for x e #(T) and _y cV (T*) 

T* exists if and only if V(T) is dense in X. The closure of P(T) in X is 
denoted by P(T), thus V(T) is dense in X implies P(T) = X. A linear operator 
T^ on V(T) C X into X will be called symmetric if T* ^ i.e., if T^ is an 

extension of T. A linear operator T^: V(X) + X is called self-adjoint if 

T = T*. 
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The following properties of symmetric and self-adjoint operators 
should be noted: 

i. A symmetric operator has a closed symmetric extension 
T** = (x*)* 5 T. 

ii. An everywhere defined symmetric operator is bounded and self-adjoint* 

iii. A self-adjoint operator is closed since an adjoint operator is 
closed. 

The relation between closed operators and their adjoints is expressed 
in the following important theorem and corollary: 

Theorem III-l . Let T be a linear operator on V (T) C X into X such that 
V (T) = X. Then tt admits a closed linear extension if and only if T** =(jC*)* 
exists, i.e., if and only if V (T*) = X. 

Corollary III-l . If t?(T) = X, then is a closed linear operator if and only 
if T = T**. 

Let us next illustrate the concept of bounded and unbounded linear 
operators on finite and infinite dimensional spaces. 

Example III-l . a. Consider the following operators in R n , 

If x = col (x^, x 2 , x^, ...» x n ) then let 

D x = col (x x , J * 2 ’ *** n X n^ 

E x = col (x. , 2 x« , ... nx). 

12 n 

n 2 l/2 

If | |x| | = ( £ xf ) then both and JJ are bounded since 

i=l 

| |D_ ac | | < | |x| | and | | jS xj | < n l |xj | . In fact, all stationary linear 
operators on a finite dimensional Euclidean space are bounded. 

b. Consider next the same operators on a infinite dimensional 
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00 

Euclidean space R , thus if x = col (x^, x 2 , • •••) then 

D x = col (x 1 , j x 2 ) 

E x = col (xp 2 x 2 » )• 

00 

If the norm of x e R°° is defined by | |x| | = ( £ then D will 

i-1 1 

still be bounded, since | |D_ x| | < | |x| | « however, JE will no longer be bounded. 

The concept of dissipative operator is defined next. 

Definition III-9 . Let T be a linear operator such that P(T) and R(T) belong 
both to the Hilbert space X. Then T^ is called dissipative if 

Re < X 2? - 0 for a H e ^ C£) • 

Definition III-10. Let T be a linear operator such that V(T) and R(T) both 
belong to X. Then T^ is called strictly dissipative if there exists a constant 
c > 0 such that 

Re < X 2£» 2? - " c I lil I ^ f° r 21 G ^(T)» 

This concept can similarly be defined in terms of semi-inner product. 
This concept is used frequently in the following chapters. Definitions III-9 
and III- 10 Imply that -T_ is positive and positive definite respectively. 

Spectral theory for an operator T^ is the distribution of the values 

of X for which T, = XI - T has an inverse and the properties of the inverse 

""A 

when it exists. The following are the definitions of the frequently used 
resolvent and spectrum. 

Definition III— 11. If X is such that the range R(T , ) is dense in X and T. 

— " — — — o — A o — A o 

has a continuous inverse (X I — T) then X is said to be in the resolvent 

o - - o t 

set p(T) of T and this inverse (X I - T) - ^ is denoted by R(X ; T) and it is 
— — o — — o — 

called the resolvent ( at X ) of T. All complex numbers X not in p (T) lorm a 

— - 1 ■ " ' o — — 

set ij(T) called the spectrum of T. The spectrum a(T) is decomposed into 
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disjoint sets P o (T), ^ 0 (T) and R o (T) with the following properties: (T) is 

the totality of complex numbers A for which does not have an inverse; P (T) 
is called the point spectrum of T. 

C (T) is the totality of complex numbers A for which T has a 

discontinuous inverse with domain dense in X; C (T) is called the continuous 

a — 

spectrum of .T. 

R^(T) is the totality of complex numbers A for which has an inverse 
whose domain is not dense in X; R 0 (T) is called the residual spectrum of 
An important theorem concerning the resolvent is: 

Theorem III-2 . Let X be a Banach space and "T a closed linear operator with 
its domain V (T) and range R(TT) both in X. Then, for any A^e p ( T) the 

resolvent (A I - T)~^ is an everywhere defined continuous linear operator. 

o — — 

The following section is devoted to some properties of Hilbert spaces, 
which are very important for the development of the stability theory for a 
class of operator differential equations. 

C. Hilbert Spaces 

In Chapter II the general operator differential equation (II-4) was 

introduced with A defined on a general normed linear space X. It is next 

assumed that X is a Hilbert space. Thus, either the norm of X is induced by 

an inner product, or the norm satisfies the parallelogram law and therefore 

induces an inner product. In Section A equivalent norms, in Banach spaces 

were defined. Since a Hilbert space belongs to a special class of Banach 

spaces, the concept of equivalent norm not only holds, but can be formulated 

more specifically with respect to the inner product structure of the Hilbert 

space. In particular, it allows the introduction of the concept c f equivalent 
inner products as defined by: 
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Definition III-12 . Let = (H, <• ,.>^) and H 2 = (H, <• **> 2 ) Be Hilbert 

spaces consisting of the elements of a linear vector space H and the inner 

products <*»»>^ and <.,.> 2 * respectively. The inner products are called 

equivalent if and only if the induced norms are equivalent. 

This concept of equivalent inner product enables one to "carry" the 

stability properties from one Hilbert space into another Hilbert space defined 

for the elements of a linear vector space. It is also possible that the norm 

induced by the inner product of a Hilbert space, H^ = (H; <*»»>^), is 

equivalent with another norm which does not satisfy the parallelogram law. 

This last norm will not induce an inner product; the resulting space (H; | | • | | 2 > 

is thus not a Hilbert space, but still a Banach space. 

The characterization of equivalent inner products can be made more 

(23) 

explicit by the Lax-Milgram Theorem . 

Theorem III-3 . (Lax-Milgram). Let H be a Hilbert space and let 'B(x»z) be a 
complex-valued functional defined on the product Hilbert space H x H which 
satisfies the conditions: 

i. Sesqui-linearity ,i.e., 

BCa^ + a 2 x 2 %) = a 1 BQ^, y) + a 2 B(x 2 , y) 
and 

B(x, 3^ + 3 2 £ 2 ) = + B 2^(— *— 2> 

ii. Boundedness, i.e., there exists a positive constant y such that 

|B(x,y)| < y ||x|| . | | Z | | . 

iii. Positivity , i.e., there exists a positive constant 6 such that 

B(x,x) > 5 ||x|| 2 

then there exists a uniquely determined bounded linear operator J3 eL(H,H) with 
a bounded linear inverse c L(H,H) such that < x,p = B(x, S. %) whenever 
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x and y. e H and | | sj | < 6~\ | |s_ _ ^| | < y . 

The implications of this theorem for the equivalent inner products are 
given by the following theorem: 

Theorem III-4 . Two inner products defined on a real linear vector space H 
are equivalent if and only if there exists a symmetric bounded positive 
definite linear operator j3 e L(H,H) such that ^c,y> 2 = < i>iL for all 
x, y e H. 

Proof: Let B(x,^) = < 2£»S. v^, where_S is symmetric. Then B(2c,x) is bilinear 

and B(£,jc) = < y_,_S x>^ = < _x,_S Since S_ is positive definite B(x,x) = 

< i>.£ > 5 < 3C and since S_ is bounded, B(jc,x) = < x.,jS _x>^ - I lil I x * 

2 

| |_S = y| lilli there follows that B(jc,;y) = < satisfies all the 

properties of an inner product. The fact that the norms induced by <*>*>^ and 
<»»*>2 are equivalent implies that the inner products are equivalent. 

Suppose <. ,.>^and <*>*> 2 are inner products with equivalent norms, 
i.e., there exist positive constants a and 3, 0 < a <8 < “, such that, 

a INI, = < IM! 2 _5 $ I M I , . 

Then since < »»«> 2 bilinear in real H, there follows 

I < 2i»Z > 2l <1 Ul l 2 l Izl 1 2 - 32 Nil I x l Izl \ V 

and therefore by Theorem III-3 there exists a uniquely determined _S such 
that 

< x,x> ± = *Q£,.S y > 2 . 

But since< = < X»i > ^ there follows < x».S y> 2 = <^,S_ 3: >2 and S_ is symmetric 

thus proving the theorem. _S is obviously a self-adjoint operator. 

In this framework of equivalent inner products it is possible to 
derive the stability properties of the general operator differential equation of 
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type (II-4) from the knowledge of the properties of the operator A only. Once 
the formal partial differential equations are formulated in terms of the 
general operator differential equations, the equivalent inner products enable 
one to derive the maximal system parameter ranges for stability, 
can be extended to complex Hilbert spaces. 


This concept 



IV. SEMI-GROUPS, GROUPS AND INFINITESIMAL GENERATORS 


In order to develop a stability theory for the solutions of the 
operator differential equation 

d x 

— = A x (x e£>(A) C X; A: V(A) -> R(A) C X) (IV-1) 

where X is a Banach space, A will be restricted to the class of operators which 
are the infinitesimal generators of bounded linear operator valued functions 
t > 0 that satisfy the condition 

T =T. *T,T 0 =I. (IV-2) 

— t+s — t —s* — 0 — 

The notion of infinitesimal generator, to be formulated later, was 

introduced by Hille and Yosida. Much about this basic concept can be found 

( 28 } ( 23 } 

in their respective books v K } , together with the general theory of linear 
operators in function spaces. Instead of the bounded linear operators 
having merely the semi-group property, the more restricted semi-groups of 
class (Cq) are introduced. It is in terms of these bounded operators that the 
properties of a system (IV-1) can be studied more easily than in terms of the 
operator A. 

In the following two sections the semi-groups of class (C Q ) , groups 
of class (Cq) and their infinitesimal generators are introduced* 

A. Semi-Groups and Groups of Class (C^) 

Let X be a Banach space with norm | |*| |. For each fixed t > 0, let 
be a bounded linear operator on X into X, eL(X,X), the set of all bounded 
linear operators mapping X into X. Then the single parameter family of 
operators {T^; t > 0 } C L(X,X) with parameter t sR + = [0, ») is said to 
satisfy the semi-group property if 
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(t,s > 0). 


Definition IV- 1. 


It 


T = T 
— s — t+s 


If {T? t ; t > 0} C L(X,X) satisfies the conditions 


group 


i- It • 2, - I t+S <'.«■ i °> 

ii, Tq = _! (I is the identify operator in L(X,X)) 

iii. lim I It x - T x| I = 0 for each t_ > 0 and all x eX, 

1 t — — t — 1 1 0 — * 

**0 ° 

then {T^.} is called a semi-group of class (C^ ) . 

(23) 

Because of the semi-group structure v 1 it follows that for a semi- 
{T^} of class (Cq) there exist constants M > 0 and 6 < 05 such that 


I ll t l I < M e 6t (t > 0) . (IV-3) 

If in addition to (i), (ii), and (iii) of Definition (IV-1) 0 in 
(IV-3) can be chosen as $ = 0, then 


III,. I I < = M (0 < t < «) (IV-4) 

and (T t > is an edui-bounded semi-group of class (Cg). In particular if 
M = 1, i.e., if 

I I— 1 1 I = 1 for (° < t < «) (IV-5) 

then {Tj.} is called contraction semi-group of class (Cq ). 

(23) 

An equi-bounded semi-group of class (Cq) is equi-continuous in t . 
The equi-continuous semi-groups are of main interest in the following- 
discussion of stability properties of unbounded differential operators. It 
is clear that since {T } is only defined for t >• 0 the corresponding operator 
differential equation does not generate a dynamical system. The system 
becomes a dynamical system if fT fc } can be extended to the negative time axis 
and {T } becomes a group: 

Definition IV-2 . If {T fc • - ~ < t< «>} C L(X,X) satisfies the conditions: 
i* i t * Is = lt+s " t » s < 
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Io = I 


iii. lim T, x - T. x = 0 for each t^eC-", + ®) and 

T “ r •“ 1 1 


t->t r 


for all x e X, then {T^} is called a group of class (C^). 

(23 ) 

The alternative definition (Yosida v ) has the group {S^} defined in 

A 

terms of equi-continnous semi-groups {T } and {T } by letting = T, f or 


t > 0 and S = T„ for t > 0, yielding S. S = I = T. T. . 

= ~-t “ t = ”t — _t “ ~ t — t 


Thus 



A 

with T^ defined on R(T ) C X. 
— t — t = 


(IV-6) 


The disadvantage of defining the group {S^} in this manner is that 
A 

as defined by (IV-6) is a semi-group only if it is defined on all of X, 


thus if R(T t ) = X. 


If (T^ is a semi-group (group) of class (C^) , then {S^.} with 

S L = e at ^ is a semi-group (group) of class (Cq) for constant a e(-«°,®). 

(23 ) 

The group structure provides the norm of T fc , where {T^} is a group 

of class (Cq) with an upper bound — there exist constants M > 0 and 6 < » 

such that 

| iTj | < M e e l C l (-«, < t < «). (IV-7) 

If in (IV-7) 6 = 0, then {T^ is an equi-bounded group of class (C^ ) . 

And if in addition M = 1 in (IV-7), then {T } is a contraction group of class 


< c o>. 


In the stability theory of semi-groups the contraction semi-groups of 
class (C Q ) are very important. The nature of the groups, however, does not 
limit the stability properties to contraction groups of class (C Q ) . An 
extremely important theorem in proving necessary and sufficient conditions for 
asymptotic stability of groups was proven by Vogt: 
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Theorem IV-1 . (Vogt). Let {T^ ; t e (-“,«>)} C L(X,X) be a group. Then there 
exist positive constants m and a such that 

I \lf. x! I > M e “I 1 "! | |x| | (-« < t < ~) (IV-8) 

Proof : Since T • T = T_ = I , 

— t - U — 


or 



* I Iz J 1 1 li, 


-t 






From (IV-7) : 



m ir 1 

> M -1 

e -3|t| 



-t 




and hence 

1 lit il 1 

5 m" 1 

e “3 | t | 

II* 


which is the required inequality. 


It should be noted that (IV-8) is a sufficient condition f or the 

existence of the inverse S ^ = S for each t e (-«,«>) with P(S ) = X. 

-t “-t -_t 


B. Infinitesimal Generators of Semi-Groups 
and Groups of Class (C^) 


In the previous section, the semi-group or group system structure has 
been established. The next step is to derive the conditions for A so that the 
solutions to (IV-1) possess this structure. In the next chapter the Lyapunov 
stability definitions and theory is established for these semi-groups or 
groups. 

Definition IV-3 . Let (T^ » t > 0} C L(X,X) be an equi-continuous semi-group 

of class (Cq). The infinitesimal generator A cf T^ is defined by 

A x = lim (h -1 (T, - I)x} (IV-9) 

h->o + 

whenever this limit exists. 
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A is a linear operator with domain 

P(A) = {x e X; lim [h -1 (T, - I)x] exists in X} (IV-10) 

~ h-0+ -- 

and the range of A, R(A) C X. Since V(A) contains at least the vector j), it 

■ ■ — (23) 

is nonempty. In fact P(A) = X 

The following theorem gives the necessary and sufficient conditions 

for A to be an infinitesimal generator of a semi-group. The proof of the 

(23) 

theorems in this section can be found in . 

Theorem IV-2 . Let A be a linear operator with P.(A) dense in X and R(A) in X 
and let the resolvent (I - n ^ A ) ^ exist in L(X,X). Then A is the infinites- 
imal generator of a uniquely determined equi-continuous semi-group if and only 
if there exists a positive constant c independent of n and m such that 

| la - n -1 A)" m | I <C ( n =l,2,3, ; m=l,2,3, . . .) . 

av-ii) 

The semi-group generated by A in the above theorem is basically 
related to the spectral properties of A in that (A_I - A) ^ exists everywhere 
for Re(A)> 0 if A generates a semigroup statisfying [ |l^_| | < M. 

Corollary IV-1 . If in Theorem IV-2, (IV-11) is replaced by 

||(I - n" 1 A) _1 || <1 (n=l,2,3, . . .) (IV-12) 

then A generates a contraction semi-group. 

The following theorem relates certain semi-groups and their 
corresponding infinitesimal generators. 

Theorem IV-3 . Let A be the infinitesimal generator of the semi-group {T^} 
satisfying 

i. T t -3. <t.« i 0). 3o - I 

ii. t lj.mQ+| jj^x “ill = 0V_xeX 
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iii. | ItJ | < M e^t (t > 0) with M > 0 and 3 < » and independent of t. 
Then (A - 3 ,1) is the infinitesimal generator of the equi -continuous semi- 

_ /5 ♦- —X 

group = e p T t cf class (C Q ) and (A I - A) is everywhere defined for 
Re[A] > (3. 

The above theorems IV-2 and IV-3 give rise to the following corollary: 
Corollary IV-2 . Let A be a closed linear operator with V(A) = X and R(A) C X 
and let the resolvent (JI - n”^ A) - ^ exist in L(X,X) for integer n sufficiently 
large. Then A is the infinitesimal generator of a semi-group satisfying (i) , 
(ii), and (iii) of Theorem IV-3 if and only if there exist constants, M > 0 and 
3 < °° such that 

| | (I - n~ 1 A)~ in | | < M(1 - n^S) - ® for m=l,2,3, ... and all large n. (IV-13) 
In particular for those semi-groups T, satisfying (i) , (ii), and 

L- 

||Tj| < e 3t f or all t > 0, (IV-14) 

(IV-13) can be replaced by 

| | (I_ - n" 1 A)"* 1 ] | < (1 - n*" 1 ^)*" 1 for all large n. (IV-15) 
Notice that 3 can be positive as well as negative in cases (IV-13) and (IV-15). 

These results relate semi-groups and their infinitesimal generators. 
Similar results hold for groups and their infinitesimal generators. The follow- 
ing theorem summarizes these results: 

Theorem IV-4 . Let A be a linear operator with V (k) = X and R(A) in X. Let 

—1 —1 

the resolvent (I -n A) exist in X. Then A is the infinitesimal generator 
of a uniquely determined equi-continuous group (T^y t e(-°° ,°°)} of class (C^) 
if and only if there exists a positive constant M >0 such that 

| | ( I - n - *A) -m | | < M (m=l,2,3, ... and all large |n|, n < 0) 

(IV-16) 
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and 


| | T fc | | = <M (-» <t<»). (IV-17) 

If in addition A is closed and there exist positive constants M > 0, g > 0 
such that 

||CX- n" 1 A)“ m ] | = <M(1 - |n _1 |B)" m (m=l,2,3, . . ,larpe|n| , n > 0) (IV-18) 
then tC ; - 00 < t < m } satisfies 

l|T t || = < Me^l (-«» <t <«). (IV-ld) 


If 


I I (I ~ n 1 A) 1 | | < (1 - | n 1 1 8) _1 (for large |n| , n > 0) 


(IV-20) 

then { - 00 < t < «>} satisfies 

| | ^ II < e 3 1 C I (— < t < «*). (IV-21) 


The statement A is a closed linear operator in Corollary IV-2 and the 

last part of Theorem IV-4 is crucial. Further elaboration is possible for the 

case when _A, defined on a Hilbert space, is the infinitesimal generator of a 

contraction serai-group of class (C^). These contraction semi-groups of class 

(C Q ) are very important for stability investigations. Their infinitesimal 

generators are the earlier defined dissipative operators. The results of the 

(24) 

following theorem and corollary are due to Phillips 

Theorem IV- 5 . (Phillips). Let A be a linear operator with domain P(A) and 
range R(A) both in the Hilbert space H and V'QO “ H, then A generates a 
contraction semi-group of class (C^) in II if and only if A is dissipative with 
respect to the inner product defined on H and R (I - A) = H. 

A consequence of this theorem is the following corollary: 

Corollary IV-3 . If A is a closed linear operator with P(A) and R(A) both in 
the Hilbert space H and l)(A) = H, then A generates a contraction semi-group of 
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class (Cq) in H if A and its adjoint A*, are both dissipative with respect to 
the inner product defined on H. 

In the corollary it is again required that A is closed, a statement not 
made in the theorem. 

The principal requirement for A to generate a contraction semi-group is 
that D(A) is dense in H and that A is maximally dissipative, i.e., A is not the 
proper restriction of any other dissipative operator. A necessary and suffici- 
ent condition is that R(X _I - A) = H for all X > 0. In this case A is closed if 
and only if R(X jt -JO is closed. However, if no conditions are imposed on 
R(X X, - A) or equivalently on R(I - A), then A must be closed. 

After having developed the relationship between the semi-group or group 
and the infinitesimal generator, there remains one more step to identify 
systems having the semi-group property with the operator differential equation 
(IV-1) . Let L) denote the time derivative of 1 x for x e X and define this 

L — L 1 1 

derivative by 

D T. x = lim [h _1 (T - T ) x] (IV-22) 

t -t - h ^Q+ “t+h -t - 

for x s X if the limit exists. 

From Yosida^^ 

Theorem IV-6 . If x e V (A) C X, then jc e V (D^T^) and 

-t - = - -t - = -t - - f° r t > 0. (IV-23) 

And in particular, A is commutative with T_ for x, e V(A). Thus if x e t?(A), 
then jC x e V(A) for t > 0. In (IV-23) of the above theorem, one easily 
recognizes the operator differential equation: 

d _x 

-r — = A x (x e V(A) C X). (IV-24) 
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Thus if in accordance with the above theorems proper restrictions are 
placed on then the solutions to (IV-24) are the semi-group or group 
trajectories given by the semi-group or group that is uniquely generated by A. 
The properties of the solutions can be studied by investigating the properties 
of this semi-group, i*e., of the infinitesimal generator A as given by (IV-24). 
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V. LYAPUNOV STABILITY THEORY FOR 


SEMI-GROUPS AND GROUPS 
A. Definitions 


Once the semi-group or group structure is established for the general 
operator differential equations 

dx 


-r- = A x (x eV (A) C X) 


(V-l) 


then the solutions are semi-group or group trajectories. Thus the solution 
starting at t=0 from ac — Xq eV(A) is given by 


x(t;xQ) = T t Xq t > 0 (V-2) 

with x(0;j{q) = x^, where A is the infinitesimal generator of the semi-group 
{T^} . If A is the infinitesimal generator of a group{ T^} , then (V-2) holds 
for -»< t< + 00 . Note that by Theorem IV-6, jc(t; Xq)e P(A) (t > 0) if Xg£t?(A). 
By the linearity of the semi-group, any trajectory can be referenced to the 
origin xfO. Thus the stability of any solution can be determined by studying 
the stability of the solution xCtjXg) = (), the origin or null solution. It is 
now possible to give the definitions of stability in terms of the semi-group 
T fc generated by A of (V-l). 

Definition V-l . The origin of (V-l) is stable in the sense of Lyapunov (with 
respect to initial perturbations) if and only if, given an e > 0, there exists 
a 6 > 0 such that 


implies that 



(V-3) 


|T t xj | < e (t > 0; V 3^ e X). 


(V-4) 
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Definition V-2 . The origin of (V-l) is asymptotically stable in the sense of 
Lyapunov (with respect to initial perturbations) if and only if 
i. it is stable 

ii. lira | | T Xq | | =0. (VXqEX). (V-5) 

t->oo 

The exponential nature of the semi-group structure usually gives rise 
to a stronger form of asymptotic stability, namely exponential asymptotic 
stability as defined by: 

Definition V-3 . The origin of (V-l) is exponentially asymptotically stable in 
the sense of Lyapunov (with respect to initial perturbations) if and only if 
i. it is asymptotically stable 
ii. there exist positive constants M and 6 such that 

I lit XqII <M e" et IIxqII (VxgtX). (V-6) 

From Definition V-3 it is clear that | |T^| | _<M e Bt . If {T^} is a 
group, then from Theorem IV-1 there follows for t > 0, | |T^| | > m e at # A 
group {T^} with the property that there exist four positive constants, 

M > 1 > m >0,a > 6 >0 such that 

m e _at | |x| | < I lit x| I ^ M e -Bt | |x| | (t > 0, x e X) (V-7) 

is called a group of exponential type . 

The following theorems are direct consequences of the above definitions 
and those in the preceding chapter: 

Theorem V-l . A sufficient condition for the stability of the null solution of 
(V-l) is that the semi-group {T^} be equi-bounded. 

Theorem V-2 . A sufficient condition for the exponential asymptotic stability 
of the null solution of (V-l) is that there exist positive constants M and 0 
such that 
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I lie I I <M e 


(V-8) 


-0t 

% 

B. Sufficient Conditions for Stability and 
Asymptotic Stability of Semi-Groups 

In the preceding section the stability properties are defined in terms 
of the semi-group properties. Thus, when this semi-group is generated by an 
operator, it is in terms of the solutions to the operator differential equation. 
However, rather than first solving the equation, one would like to base the 
stability properties of the system directly on the properties of the operator, 
i.e. , to determine conditions for the operator so that the solutions exist and 
at the same time are stable. These conditions are clearly spelled out in the 
theorems in Chapter IV, Section B. 

It is intuitive that stability of a null solution requires A to be the 
infinitesimal generator of an equibounded semi-group. However, it is not so 
easy to relate asymptotic or exponential asymptotic stability directly to such 
a basic property. However, once the contraction property is established (IV-14) 
and (IV-15) of Corollary IV-2 seem to provide the answer. 

In Theorem IV-5 the principal conditions on A to generate a contraction 
semi-group are that A is a dissipative operator with respect to an inner 
product K *»* > ^ and R(_I - A) = H^ and consequently the null solution is stable. 

The dissipativity is defined with respect to the particular inner 
product of the space, but the semi-group property is invariant under equivalent 
norming. In general, a semi-group { T^ is stable if it is equi-bounded, i.e., 

| lij. ||i _< The invariance under equivalent norming does not necessarily mean 
that if | |T t | ^ < M 4^ | l^l I 2 = 1 then I 1*1 1 2 corresponds to an inner product 

<•*•> 2 . In other words, | |.| tpay define a Banach space- rather than a Hilbert 
space. 


47 



The invariance of stability under equivalent norming suggests that if 
A is dissipative with respect to any inner product equivalent to the inner 
product of the space, then the dissipative property of A with respect to this 
inner product is sufficient for stability and generation of a contraction 
semi-group. For this reason the dissipativity of A is extended as follows: 
Definition V-4 . Let be a real Hilbert space with inner product <»»«>^. 

A is dissipative in the extended sense if and only if there exists a self- 
adjoint, positive definite, bounded linear operator S_ eL(H,H) such that 

< x, _S A x>! < 0 x e P(A) (V-9) 

and A is strictly dissipative in the extended sense if and only if there exists 
a g > 0 such that 

< 2 £> A. x>^ < -g < _x » S_ x>^ x zO (A) (V-10) 

Since under the conditions on S_: 

< x, i >2 = < x.» S. v> x x, 2 e H 

one could also call A dissipative in the extended sense if there exists an 
equivalent inner product with respect to which A is dissipative. The following 
theorem follows directly from Theorem IV-5. 

Theorem V-3 . Let A be a linear operator with V(A) = H and R(A) C II and 
R(I - A) = H. Then A generates a contraction semi-group if and only if A is 
dissipative in the extended sense. 

Proof : Since A is dissipative in the extended sense there exists an inner 

product with respect to which A is dissipative. The "if" part follows from 
Theorem IV-5 by using this inner product as the inner product for the space. 

If A generates a contraction semi-group according to | |»| | 2 = 
then A is dissipative with respect to <•»•> and hence is dissipative in the 
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extended sense. 


Remark : Theorem V-3 includes Theorem IV-5, in the sense that A is not 

dissipative with respect to the original inner product of the space. In' this 
case {T^; t > 0} is at least equi-bounded with respect to the norm induced by 
the original inner product. 

In order to derive some principal results on the stability and 
asymptotic stability of certain classes of semi-groups, Lyapunov functionals 
will be introduced. 

A Lyapunov functional on a real Hilbert space is defined through the 
symmetric, bilinear form, 

V(x*1L) = < x* S. JL x.» X e H (V-ll) 

where S is a self-adjoint, bounded positive definite, linear transformation, 
jS eL(H,H). Define the Lyapunov Functional by 

v(x) = V(x,x) xe H. (V-12) 

The time derivative of v(x), denoted by v(x) along solutions to (V-l) 
with A generating a semi-group {T^ ;t > 0} C L(H,H) is given by: 

v(x) = lim ~ (V(T x, T x) " V(x,x)) (V-13) 

t->0 c z 

whenever this limit exists. 

But with 

V (T t x + x, T fc x - ^c) = V (T fc x, T t x) + v (x> ~ v <£ t x >x) ” v (x» x) 

and since V(x,%) is symmetric it follows that: 

v(x) = lim y V ( (T + I) x, (T - I) x) 
t-*0 

= 2 V (x, A x) (x s 5(A)). (V-14) 

This leads to an important result: 
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Theorem V-4 . Let A be a linear operator with $ (A) = H, R(A) C H, and 

R(JC - A) = H. Then A^ is an infinitesimal generator of a contraction semi-group 

if and only if there exists a Lyapunov Functional v (x) such that 

v (x) ~ 2V (x, A x) < 0 x e V (A) . 

Corollary V-4.1. Under the conditions of Theorem V-4, the null solution of 
(V-l) is stable. 

Proof : The "only if" part cf the theorem follows from Theorem V-3 by taking 

v(x) as, 

v (jc) = < 3C,X>. 

If there exists a v(x) given by V(ic,^), then let < x»2L > 2 = V(jc,£), with 
< x> A x >2 <0. Thus by Theorem V-3, since V(A) = H and R(I_ - A) = H, _A 
generates a contraction semi-group t > 0} with 

Mltll 2 <i t > o. 

This result implies Corollary V-4.1. 

The statement of the theorem indicates that the conditions for existence 
and stability of the solutions to a general operator differential equation of 
the form (V-l) are much more restrictive than when (V-l) represents a finite 
dimensional system of ordinary differential equations. The free interchange 
of equivalent inner products can facilitate the investigation of these require- 
ments. 

The following theorem points out its importance for the case of 
asymptotic stability. 

Theorem V-5 . Let A be a linear operator with V(A) = H, R(A) C H, R(I^ - A) = H. 

Then the null solution of (V-l) is asymptotically stable if there exists 

Lyapunov Functional v(x) such that 

v(x) = 2V(x, Ax) £ - Y | | x | 1 1 x e ^(A). 
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Proof : Since all the conditions cf Theorem V-4 are satisfied it follows that 

A generates a contraction semi-group. To show that the null solution is 
asymptotically stable notice that from 


certainly follows 


v(x) = 2V(x, A x) < - y | |x| \ \ 


V (x, A x) = < x , A x>„ < - a | | x | L . 


Next it must be shown that for n > H > 0, the inverse of (I - n~ A) exists such 


| | (I - n" 1 A)" 1 | | < (1 - n “ 1 e) _1 


f or some g. From the relation 


< (_I - n ^A) = < ijiSr’o ~ n ^ < A x, x>^ 


it follows from Schwarz's inequality that 


(I, - n _1 A)x| | 2 I |x| l 2 > I <(I - n -1 A)x,x> 2 1 = 


i < A x,x> 0 

i •. —1 — z 

= < x,x> 0 1 - n j 


and therefore that 


It follows that 


| |(I- n _1 A) _1 | L < (1 - n” 1 (- a))" 1 n > a. 


m —at 

l-t I 1 2 = e 


This is a sufficient condition for the asymptotic stability of the null solu- 
tion. 

Remark . The above theorems can similarly be stated with the condition of A 
being closed and with the subsequent modifications following the theorems of 
Chapter IV. 
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Note also that the Lyapunov Functional gives only sufficient conditions 
for stability or asymptotic stability respectively. This is a distinct 
difference between the stability properties of systems having the semi-group 
and group properties respectively, a fact clearly demonstrated in the follow- 
ing section. It is foreseen that a slight conceptual change in the selection 
of the Lyapunov Functional may alleviate this difficulty. 

C. Necessary and Sufficient Conditions for the Exponential 
Asymptotic Stability of Groups 

In the preceding section the sufficient condition for the asymptotic 
stability of semi-groups in terms of its infinitesimal generator was establish- 
ed. This applies also to groups. However, the group property is the funda- 
mental property of dynamical systems for which a Lyapunov stability theory has 
been developed which includes both sufficient and necessary conditions for 
stability or asymptotic stability. Thus, it is not surprising that in the case 
of groups the approach developed in the previous section can be extended to 
include the necessity. Again the Lyapunov Functional not only gives necessary 
and sufficient conditions for stability or asymptotic stability but guarantees 
also the existence of the solutions. 

Theorem V-6 . Let A be a linear operator with P(A) = H with inner product 
and such that I - A) * H for real X and|x| sufficiently large. 

Then A is the infinitesimal generator of a group of exponential type if and 
only if there exists a Lyapunov Functional v(x) = V(x,jO with 

V(jc,^r) = < 3c, _S y>^ = < jc, where S is a symmetric or self-adjoint bounded, 

positive definite operator, J3 e L(H,H) and such that for some constants a and 
3, a > 3 > 0 
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- 2 <xV(x,x) <v(x) = 2V(x, Ax) < -2 B V(x,x). (xe£>(A)) (V-15) 


Proof : If there exists a V(x,£> with the above properties, then let 

<•,.>2 = V(.,.). To show that A generates a group there follows: 



< ( I_ - n 1 A) 2 £»x > 2 

= < 

2£>x > 2 ~ 

n 1 < A x, x> 2 • 


Then with 

Schwarz's inequality: 






| | (I. - n ^xl 

ll 2 i 

1 1*1 1 2 

, < A x,x> 0 
| i - n -1 "- 2 -| . 

< x»2i > 2 


For 

n > N g > 0, 

III-”' 1 a|| 2 > 

1 + 

I”' 1 ! 6 


(V-16) 

and for 

n < - N 
= a 

< 0 





III - n_1 A|| 2 > 1 - In'* 1 ! a. 

Thus for |n| > max (N^, N^) n < 0 

I \L “ n -1 A | | 2 > 1 — In” 1 | a > 0. 

Therefore (Jt - n 1 A) 1 exists and 

| | (I - n" 1 A)" 1 1 1 2 < (1 - |n _1 |ct) -1 

and thus by Theorem IV-4 A generates a group {T^_; t e (-“, “) } which satisfies 

| | Tj. | 1 2 < e a I Z I (-00 < t < «) . 

From inequality (V-16) it follows by Corollary IV-2 that A^ generates a semi- 
group {T^; t > 0} with 

I lltl I2 = e ~ Pt t e t°» “>• 

Since ||«||^ and ||.|| 2 are equivalent, this completes the proof of suffi- 
ciency, The only if part will be demonstrated by construction of the required 
Lyapunov Functional, If A generates an exponentially asymptotically stable 
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group, then it follows from Theorem IV-1 that the most general group 
{T^; te(-°°,«>)} satisfies 

® e ~ at | lit \i * I llj. il \i * M e_Bt I lilli te[0,oo) 

where »>M>l>m>0, «»>a>8>0 and | |t | < M e 0 !^! (-» < t < »). 

Next take n 

V(n; x,2.) » / < l t x» It £>1 dt (V-17) 

with n >0. Since ^T^x* I t X > 1 ds f° r eacdl t a numerical value, which is 

everywhere defined and continuous on the compact interval [0,n] (V-17) can 

be interpreted as a Riemann integral. For each fixed n, V(n; x»X) *- s a 

symmetric bilinear form, which satisfies 

|V(n; x»X) I < I |x| | x I Ixl l x M 2 / e ~ 2 ^ dt 

o 

- Mil l x llxlli fs (1 - e " 2en) - ( v - 18 > 

Next let V(x_,x) - lim V(n; then from (V-18) it follows that V(jc,£) 

n-x” 

exists and is well defined, moreover 

l v Qi»X) I < fg I |x| li Mil li 

2 

and V(x,x) > —■ | |x| \ 1 . 

The symmetric bilinear functional V(x,y) satisfies all the conditions 
of the Lax-Milgram Theorem, thus there exists a uniquely determined, bounded, 
Symmetric, positive definite linear operator with 

V(x,£) = <x, j> 2> 1 ~ < £ x, ^ > 1 x, Z eH 


with 


M 


llslli < 28 * 

Next let v(x) = V(jc,x). Then there remains to be proven that v (x) 


V(x, Ax) 


satisfies the conditions of Theorem V-6. But 
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V&t- 2k* It x) - V(x,x) 


lim 

n->» 


/ < T T. x, T T x>- dr 
i — r— t — * — r — t — 1 


n 

- / < Tx, T x>- 
k — t— — r — 1 


dx 


Therefore 


= lim 
n->«° 


= lim 
n-*» 


n+t 


/ < T x, T x>- ds - / < T x, T x>, ds 

{_ — s — ’ — s — 1 £ — s — ’ — s — 1 

t 

- f < T x, T x>, ds. 
; — s— — s — 1 


n+t 

f < T x, T x>, ds 
n 


2V(x,A x) = lim ■— 

t->0 


VCIt 2£» l t x> " V(x,x) 


x t 

= - lim — J < T x, T x>. 

t L -s -s - 1 


ds 


t->0 


= - < x»2L>i 


(x e P(A)) 


But since is equivalent to = V(.,.), it follows that there exists 

a > 3 > 0 such that 

- aV (x,x) < V (x, Ax) < - 3 V (x.,x) . 

This completes the proof of the theorem. 

The interesting aspects of this theorem become clear in Section VIII, 
where certain wave equations are studied. There, the difference between the 
stability properties for the semi-group and group structures are demonstrated 
by taking the same simple operator studied in Section VII in a diffusion 
equation setting, and in Section VIII in a wave equation setting. In the 
next section, the formal partial differential operator is synthesized to give 
the appropriate representation (V-l). 
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VI. FORMULATION OF FORMAL PARTIAL DIFFERENTIAL OPERATORS AS 


UNBOUNDED OPERATORS AND THE SELECTION OF NORMS 

In this chapter an unbounded operator A is obtained from a class of 
formal partial differential operators. This means that V (A) and:R(A) are 
determined along with the transformation A x = ^ for certain x« 

This is carried out by means of a synthesis process. First the linear 
space in which the operator is defined is completed by introducing generalized 
derivatives. Some basic definitions and lemmas from the theory of distribu- 
tions are given. Sobolev spaces are introduced; these spaces are significant 
for the determination of the domain and range of the unbounded operator. The 
Sobolev Imbedding Theorem is stated. The chapter concludes with a discussion 
of some aspects of the selection of norms for the stability problem. 

A. Derivation of the Differential Operator from the 
Formal Partial Differential Operator 

The derivation of the formal partial differential operator in the 

(29) 

proper space setting is adapted from Dunford and Schwartz . The real case 
is considered only; the results can be extended to the complex case. 

1. Notations 

The symbol J will denote an index , i.e., a k-tuple J = (j^» • • • » j^) of 

integers, |j| = k, min J = min j., max J = max j.. The case J is 

1 < i < k 1 1 < i < k 1 

vacuous is denoted by J =0. The symbol R n is the real Euclidean n-space. 

An index J will be said to be an index for R n if min J > 1 and max J < n. If 
jc eR n , so that x. = (x^,X 2 , . .., x^) and J is an index for R n , so that 
J = (j^, ..., j^), k = |j| , then ^ will denote the expression xj^ xjg ••• x-j^. 


56 



3 3 

The operations * and r— of partial differentiation are sometimes 

o d S 

written as 3 x or 3 and 9 g respectively. If J is an index for R n and |j| = k, 
then the higher partial derivative 



9 x . 9 x . . . ,9 x . 

3 1 3 2 J k 


will be called a partial derivative of order k = |j| and will be written 9 J . 
If | J | =0 the operator 9^ is defined as the identity operator. 


2. Formal Partial Differential Operator 

If m is a positive integer, an expression 

t = 7 a. (x) 9 J 

| j| <m J “ 

where the real coefficients aj are infinitely differentiable functions in an 
open set & C R n , will be called a formal partial differential operator defined 
in ft and m will be called the order of t. 


x*(0 = l (~1) J 9 J [a (x) (•)] 

| j| < m J ~ 

is called the formal adjoint of t. In particular, if x = x* then x is called 
formally self-adjoint . If 


A 

X 




J 


is another formal partial differential operator defined in ft, then for func- 
tions f which are infinitely often differentiable in ft, x (xf) is also defined 
in ft. This is written 


x x = x = £ b'y (x) 9 J . 

|j|< m+m 
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Similarly the sum of x and x is defined by 

x = T + T= £ (a (x) + a (x)) 3 J . 

| J| < max(m,m) 

Example VI-1 . Let x. = (x-^x^x^), then a particular t can have the following 
form: 


T = V 


U\< 


<4 


ij(x) 3 J = a.^ (3c) 3 ^ + a^te) 3 ^ + a 2 (x) 3 ^ + a^(x) 3 + a Q^ 3 


' h ' 73 


+ a. 


<x) 


n?l — 2 

3x5x„3 x„ u> ’ a x„ 3 x. 


ju ^ 2 u 


2 3 


3 


3 2 

+ a 2 o 0 ^ 2 ~ a l 1 0 ^ 

» * 3xj * ’ 3 x 1 9x 2 


+ a Q (x) . 


Notice that the elements a^Cx), 0 <|j|< 4 can be zero, as is a^te). In most 
practical applications, some components of a^Cjc), 0 <|j|,< 4, will be zero. 


t*(-) = l (-1) J 3 J [a (x)(-)] = a ~r [a n (x)(0] - 


j| = <4 


3 x 


4 1 0,4,0"-' 


[a. - . (x) (•) ] - 


. 2 , . L 2,1,1"-' 

3 x^ 3 x 2 3 x^ 


2. [a 0,2,l ( - )( * )]+ 


3 X 2 3 X 3 


+ 



t a 2,o,o«<^ 


3j 

3 x 1 3 x 2 


[ a l 1 Q (x) (•) ] + a Q (x) . 


3. Function Spaces C^(fi) and Cg(ft) 

Let ft be an open set in R n and ft its closure. Then the set of scalar 

functions f defined on ft with all partial derivatives of order not more than 

k k 

k existing and continuous is denoted by C (ft) . The set (ft) consists of 

1c ^ w 

those functions in C (ft) which vanish outside a compact subset ft* =» ft’, ft’ is 
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1C 1c t_ ™ 

a proper subset of ft. Thus C^(ft) C C (ft). The set C K (ft) consists of all 
functions defined on ft having all partial derivatives of orders up to k inclu- 
sive at each point of ft and such that each partial derivative has a continuous 
extension to ft. If this is the case, 9^f(3c) is defined for x e ft and |j| < k 

T — ^ 

as the extension by continuity of 3 f(x) from ft to ft. We accept CQ(ft)=C Q (ft); 
if f e C k (ft) = C k (ft), then 3 J f « 0 on ft - ft. 

oo 00 

The spaces C (ft) and Cq ( ft) are especially important; 

4. Norm in C k (ft), C k (ft) 

k 

C (ft) is made into a F-space (a space which does not satisfy 
I I “ill I = |°| | |x| I °f the previously defined norm properties; I | • | | ^ so defined 
is called a F-norm) as follows. Let be an increasing sequence of compact 
subsets of ft or ft. Suppose that is such that any compact subset of ft 
belongs to one of the sets K^. Then for a function f in one of the spaces 
C k (ft), C k (ft), place 


p(f; J,m) = sup | 3 J f(x)| 

x eK 

— — m 

and define the F-norm of f by 

ll*llc k " l I ,1, ih- - ^ J, g ) - < VI -w 

m=l j-0 |j|=j 2 m 2 3 j! 1 + y(f;J,m) 


This norm makes the space into a complete F-space. If k < 00 and ft is compact, 

k - 

but not otherwise, the space C (ft) is a Banach space under a norm equivalent 

to the F-norm such as I 1. 1 Lk,^.. = sup I 3^* I . It is in the sense of these 

C W | J| <k 

1c 

norms that reference is made to the topology of C (ft) . 

5. Distributions 

It is essential to apply partial differential operators to complete 
spaces, i.e., to spaces with elements which all possess the necessary 
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differentiability requirements in order to carry out all required mathematical 
operations properly. Suppose, that the operator 

(t o £) - r ^ T x ^ <s) £ e r2 

, 00 2 

defined for each function f in Cq (R ; is considered. This operator Tq is 

2 2 2 
densely defined in L (R ) , the space of all functions f such that | f Cx ) | is 

Lebesgue integrable over R n . But is not closed. Let t be its closure, 

then P(t) contains nondiff erentiable functions. Which non-diff erentiable 

functions? One might expect the answer to be those (non-diff erentiable) 

2 2 

functions f such that 3 3 f belongs to L (R ). Thus one should be able to 

X 1 x 2 

define 9 9 for every function, differentiable or not, and irrespective 

X 1 x 2 

2 2 

of whether 9 9 f belongs to L (R ) or not. Such a "derivative" can no 

X 1 x 2 

longer be an element of any space of functions, but can only be a "function" 
in some generalized sense. Laurent Schwartz has developed a theory of such 
"generalized functions" in his theory of distributions. 

Definition VI-1 , i. Let ft be an open set in R n . Let {<t> n ) be a sequence of 

00 OO 

functions in Cq (ft) and let <f> be in Cq (ft) . Then if there is a compact subset K 
of ft such that all functions $ vanish outside K, and if in addition, <J> n + <J> in the 

OO 

topology of Cp(ft) we shall write 

$ n x $ in 

ii. A linear functional F defined on c”(ft) such that F(<}> n )-»-F(<f>) 

whenever 4> -»■ <j> in ft is called a distribution in ft. 

n " " 

iii. The family of all distributions in ft is denoted by D(ft). 

Definition VI-2 . Let ft be an open set in R n . Let f be a function defined in 
ft which is (Lebesgue) integrable over every compact subset of ft, then the 
distribution F, defined by 
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F(<|>) = f<p (x) f (x) d x <i> e c“(fi) 
ft 

is called the distribution corresponding to f . 

If a distribution corresponds to two functions f and g, then f = g 
almost everywhere. 

Definition VI-3 . A distribution F which corresponds to a function f in the 
sense of Definition VI-2 is said to be a function. If f is continuous, 

O Tr OO 

differentiable, belongs to L (ft), C (ft), Cg(ft), etc., F will be said to be 

2 k 00 

continuous, differentiable, belong to L (ft), C (ft), C Q (ft), etc. respectively. 

Thus a distribution, which is a functional, is identified with the 
function to which it corresponds. The introduction of distributions enables 
one, as will be seen later, to formulate a partial differential operator in 
a complete function space, in particular a Hilbert space. The following 
definition shows how a distribution may be differentiated partially. 

Definition VI-4 . Let t be a formal partial differential operator defined in an 
open subset ft of R , and with real coefficients in C (ft) . Let F be a distri- 
bution in ft. Then x F will denote the distribution defined by the equation 

(xF)(<f>) = F(t*<(>) 4> eC~(ft). 

The fact that ( J> n ^ < { ) implies x*<j> n ^ x*<j> is the justification for this definition. 
Additional justification for Definition VI-4 is provided by the statements of 
the following lemma: 

Lemma VI-1 . Let ft be a subset of R n . 

i. If the distribution F in ft corresponds to the function f in C n (ft) 
and if x is a formal partial differential operator of order at most n defined 
in ft, then xF corresponds to xf. 
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ii. 

t(cxF + 6G) = axF + 6 tG 


F, G eD(ft) 

iii. 

(a + Bt 2 ) F = a ( t jF) 

+ 6(t 2 F) 

F eD(ft) 

iv. 

( T 1 T 2 )F = t 1 (t 2 F) 


FeD(ft). 


This lemma provides the -justification for addition, multiplication by 

a scalar, differentiation, etc., of distributions. Additional properties can 

. - . . (29) 

be found m 

6. The Sobolev Space H^(ft) 

The Sobolev space H (ft) constitutes an important class of subspaces of 
I) (ft), which are at the same time spaces of functions. The following lemma is 
elementary : 

Lemma VI-2 . Let ft be an open subset of R n and let F be in D(ft). Then F is 
2 

a function in L (ft) if and only if there exists a finite constant K such that 

|F(*)| < K| |<j>| | 2 , 4. eC~(ft) 

2 2 

where | | <j> | L = /<t> d x with the integral taken in the Lebesgue sense. 

1 ft 

Definition VI-5 . Let ft be an open subset of R n and let k be a non-negative 
integer . Then 

J 2 

i. the set of all F in D(ft) such that 3 F is in L (ft) for all 

ii k k 

| J | < k will be denoted by H (ft) . For each real pair F,G in H (ft) we write 

<F » G> k = J, / 0 J F)(x)(3 J G)(x)dx 

K | J| <k ft 

and 

l|F|l k - (< F,F> k ) 1/2 . (VI-2) 

1c 1c 

ii. the symbol H^(ft) will denote the closure in the norm of H (ft) of 

OO 

the Cg(ft) functions. 
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II 

Lemina VI-3 » Let ft be an open set in R * Then the space H (fi) of the 

1c 

preceding definitions is a complete Hilbert space, and the space Hg(ft) is a 
closed subspace of H (ft). Moreover: 

Hg(ft) = H°(ft) = L 2 (ft) 

H k+1 (ft) C H k (ft) k > 0 

lrj.1 If 

H 0 C Hj(fi) k > 0. 

Lemma VI-4 . Let x be a formal partial differential operator of order k and 

let j > k. Let ft^ be a proper subset of ft, an open set in R n , such that ft^ is 

00 

compact and let a (x) belong to C (ft) , then F+tF where FeH J (ft Q ) implies 
J ^ 

1-k 

xFe H J (ftp) and the mapping is continuous. 

F e HpCftp) implies xF e Hq (ft^) and the mapping is continuous. 

7. Sobolev Imbedding Theorem 

Definition VI-6 . Let p be a point of the subset A of R n . Then A is said to 

be smooth in the vicinity of p if there exists a neighborhood U of P and a 

mapping <J> of U on a spherical neighborhood V of the origin such that 

i. (J) is one-to-one, <J> is infinitely often differentiable, and <j> - ^ is 
infinitely often differentiable. 

ii. <J>(AV) =vnixe R n |x^ = 0}. 

If the set A is smooth in the vicinity of each of its points, it is 
said to be smooth , or to be a smooth surface . 

Theorem VI- 1 . Let ft be a bounded set in R n . Suppose that is a smooth 

surface and that no point of the boundary of ft is interior to ft. If 2k > n 

and 0 < m < k - ^ then every derivative of order not more than m of the 

k 1c 

distribution F e H (ft) is continuous and the imbedding operator from II (ft) 

into C m (ft) is bounded and completely continuous. 
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8. Differential Operator Representation 

In order to apply the stability theorems of Chapter V, the formal 
partial differential operator must be extended to he defined in a Hilbert space, 
i.e., in a complete normed linear space, with the norm induced by an inner 
product. The important requirement is that the space of differentiable func- 
tions is complete, i.e., that the operator, its domain and its range are 
defined in such a way that all required formal mathematical operations are 
valid within the setting of the space. 

Let t be a formal partial differential operator defined in a domain ft 

2 

of R n , then Tq(x) and T^(x) will denote the operators in L (f2) defined by the 
equations 

t?(T 0 (x)) = c“(ft) ; T 0 (x)f - Tf f eP(T 0 (T)). 

2?(T 1 (t) ) = {f e D(fl)| feL 2 (ft), xfeL 2 (ft)}; (VI-3) 

T^x) f = x f f e DCT-^x)). 

Then by Definition VI-1 there follows that T^(x) = (T^(x*))*. Thus by Corollary 

2 

III-l, T. (x) is a closed operator in L (ft) . Moreover by Lemma VI-1 T.(x)C T, (x) 

JL 0 13 1 

is a closed extension of T^(x). 

2 

Thus T^(x) operating on an element of ^(T^x)) C L (ft) will result in 
2 2 

an element in L (ft). L (ft) is a complete linear space. The proper Hilbert 
space setting is obtained by defining a norm, [|.|[ 2 “ ||»|L, on L 2 (ft) based 
on the inner product: 

< f »g> n = / f g d x f, g e L 2 (ft) . 

U ft 

and ||f|l 0 - (< f,f> Q ) 1/2 . (VI-4) 

In order to derive the unbounded operator A from T^(t) a domain and 

2 

range in the Hilbert space L (ft) must be specified. The requirement of the 
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existence of bounded derivatives and tbe necessity for a proper adjoint 
relationship for A in the Hilbert space setting dictates a linear subset of 

O 

L (ft). These requirements can be met by specifying the domain as a Sobolev 
space. 

For a formal partial differential operator, this could lead to a 

2 

specification for A as: the operator in the Hilbert space L (ft) defined by: 

P(A(t)) = 0(A) = V(T z (r ))f) Hg(ft)H H 2k (ft) 

A f = T^t) f , f e P(A). 

2 

Under these conditions R(A) will also be in L (ft). 

Another way of achieving the adjoint relationship for A is by defining 

generalized derivatives in the sense that integration by parts formally carries 

through^ 2 ^. In both cases it places restrictions on the values of the 

distributional derivatives or derivatives on the boundary. Thus if ft=ft +3 ft 

lc 

has a sufficiently smooth boundary, 3ft, then the space H^(ft) has (distribu- 
tional) derivatives up to order k-1 which approach zero as the boundary is 
approached in ft. 

A more detailed discussion concerning this formulation of differential 

( 29 ) 

operators in the appropriate spaces can be found in Dunford and Schwartz . 
This development allows us to apply the stability theory developed for operator 
differential equations in Chapter V to a class of partial differential equations 
that can be placed in the above Rilbert space setting. 

B. Selection of Norms 

i 

It has been pointed out that stability is defined with respect to a 
norm. However, it is obvious that in formulating the particular stability 
problem the norm cannot be chosen at random. In many instances one is able 
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to select a norm in a natural way by considering the physical properties of the 

system. This leads for the wave equation to a consideration of the energy 

2 

which provides as the natural norm the L -norm as defined by (VI-4) i.e., a 
Hilbert space setting. 

One the other hand, for the heat equation the natural norm is given by 
the supremum of the temperature, thus (VI-1) with J = 0. Other examples can 
be given where the natural norm is not induced by an inner product or does not 
satisfy the parallelogram law. 

The stability theory is based on the Hilbert space theory of dissipa- 
tive o perators , the generators of contraction semi-groups or groups . Thus the 
Hilbert space structure is essential for the stability theory. 

The specific problem considered here concerns stability with respect 
to Hilbert space norm. Whenever the natural norm does not correspond to an 
inner product one can only draw conclusions about the stability properties 
with respect to this norm if it is equivalent with the particular Hilbert space 
norm for which the stability is determined. 

The role Sobolev’s Imbedding Theorems might play in relating stability 
properties under different norms, not necessarily equivalent, is suggested for 
future research. 
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VII. STABILITY OF AN EVOLUTION EQUATION 


In the previous chapter it was pointed out that in order to define A 
as an operator on a complete space, it is necessary to introduce distributions. 
A_ is then defined on some Sobolev space. A is thus associated with generalized 
functions and the conditions for stability of the solutions of (V-l) must be 
interpreted as stability conditions for the distribution solutions. 

The Sobolev Imbedding Theorem makes it possible to determine- if and 
when the stability properties can be extended to the -norms. 

The Dirichlet boundary value problem for strongly elliptic partial 
differential operators will be studied to demonstrate this -development. Section 
A is devoted to some basic properties of strongly elliptic formal partial 
differential operators. The basic theorems concerning the formulation of the 
Dirichlet problem in the distribution sense conclude this section. 

The actual stability problem is solved in Section B. 

A. Elliptic Partial Differential Operator 
Definition VII-1. Let 

t = l a (x) 3 J 

m<p 

be a formal partial differential operator of order p with real coefficients 
defined in a domain 0 of R n . Then if for each nonzero vector e R n we have 

1 a _ (x) 5. J $ 0 x e ft 

|J|=P 

the operator t is said to be elliptic . 

The boundary value problem to be discussed considers the strongly 
elliptic partial differential operators, i.e., those operators of even order 
for which Garding's Inequality holds. 
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Lemma VII-1. Let t be an elliptic operator of even order 2p with real 


coefficients defined in a domain Kq of R n . Let SI be a bounded open set such 
that SI C fig. Let 

T= l a (x) 3 J (VII-1) 

|J|<2P J_ 

and suppose that 

<-l) P l a (x) i 1 > 0, It S ( 11, n R n . (VII— 2) 

[J|-*P J " “ “ 0 - - - 

Then there exist constants K < 00 and k > 0 such that 

< xf,f> 0 + K< f,f> 0 > k | |f||2 f e c“(n), (VII-3) 

or <(t + KI ) f,f>Q > k ||f|| p . For K=0 we get 

< Tf,f> 0 > k| | f | 1^ = k< f,f> p 

or < -xf,f>Q < - k < f,f> p < - k< f,f>Q. 

x is called strongly elliptic and (VII-3) is known as Garding's 
Inequality. 

The next step is to specify a differential operator for x (still without 
boundary conditions) which later might be a candidate for the closed extension 
of an operator associated with the Dirichlet problem. T^(x) is the closed 
operator defined by (VI-3) . 

Theorem VII-1 . Let x be an elliptic formal partial differential operator of 

even order 2p, satisfying the hypothesis of Lemma VII-1, and defined in a 

bounded domain ft in R n which satisfies the hypothesis of the lemma. Let 

2 

T = T(x) be the operator in the Hilbert space L (S2) defined by the equations 

P(T(x)) = t?(T) = P(T 1 (t))0hP («)) (VII-4) 

Tf = T 1 (x)f, f e P(T). 
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Then a(T), the spectrum of T, is a countable discrete set of points with no 
finite limit points and for X i o (T), R(X; T) is a compact operator. 

The restriction of the domain of T to a subspace of the domain of 
T^(x) by (VII-4) makes R(X; T) a compact operator. 

Corollary VII-1 . Let the hypotheses of Theorem VII-1 be satisfied. Then there 
exists a constant K < 00 and a constant k > 0 such that 

< Tf,f> 0 + K < f,f> Q > k | |f | , f e Vi?) (VII-5) 
The next step is to specify the adjoint. 

Theorem VII-2 . Let x be an elliptic formal partial differential operator of 
even order 2p, defined in a bounded domain fi. Suppose the hypotheses of Lemma 
VII-1 are satisfied. Let T and S be operators in the Hilbert space L2(ft) 
defined by 

vi?) = vi? ± iT))n hJJgi), v(s) = vw^nn h£ w 

Tf = T 1 (x)f, f eV(T); Sf = T^x^f, f eVi S). 

Then T = S* and S = T*. 

Moreover we have thus < Tf, g>Q = < f, Sg>Q for any f e P(T) and 
g e P(S). 

The Dirichlet problem can be formulated as follows: 

Definition VII-2 . Let SI be a domain in R n whose boundary contains a part 

r which is a smooth surface. Suppose that no point in r is interior to . 

k— 1 — J 

Let k be a positive integer. Then if f is in C (fi) and 3 f (x) vanishes for 
all x e T and all J with |j| < k - 1 we will say that f satisfies the 
Dirichlet condition of order k on T or that f and its first k-1 normal 
derivatives vanish on r and we write: 
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(VII-6) 


O v (D) j f(x) = 0 X e r, 0 < j < k-1 

Remark ; The subscript v in the preceding formula indicates the derivative in 
the direction of the normal to r. If r is a closed rectangular (n-l)-dimen- 
sional hypersurface in R° with side I\ perpendicular to the x^-coordinate axis 
then (VII-6) becomes 

9 j f (x) 0 < j < k-1 

t — =0 x e 

3 x~ 1 < i < n. 


The following theorems will identify the operator x with Dirichlet 

boundary conditions with an appropriate operator A defined on distributions and 

thus closed. The notation T(x) is used for the operator in the Hilbert space 
2 

L (ft) defined by the equations 


£>(T(t)) = Hjj(ft)0 H 2p (ft) 


(VII-7) 


T(x)f = if ; f e P(T(x)). 

From Theorem VII-2 and the fact that Hg(ft)flH^(ft) is a dense subset 

2 

of L (ft), it follows that T(x) is a closed operator. 

Theorem VII-3 . Let x be a strongly elliptic formal partial differential opera- 
tor as given by (VII-1) and satisfying (VII-2) . Let ft be a bounded subdomain 

ft C ftp. Let 9 ft be a smooth surface and let no point in 3 ft be interior to the 

A 2 

closure of ft. Let T and T be operators in the Hilbert space L (ft) defined by 
the equations; 

P(T) = V(T) = {f e C”(ft) | f(x) = 

(VII-8) 

= 3 v ( 3 ft) f (x) = ... = 9 P“ 1 ( 3 ft) f (x) = 0,x e 3 ft} 

A A 

Tf = xf , Tf = x*f, f e P(T) = P(T). 
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Then let W and W denote the operators whose graphs are the closures of the 

A 

graphs of T and T respectively. Then 
A A 

i. W* = W , W* = W. 

ii. o( W) is a countable discrete set of points with no finite limit 
point. 

iii. It \ i o(W), R(A; W) is a compact operator, 

iv. If \ i a(W), R(A; W) is a continuous mapping on rf^ft) into 
rf n+2p (JJ) for every m > 0 . 

v. If Wf is in H m (ft) and if m > - p, then f is in C p ^(ft) and f 

satisfies the boundary conditions defining P(T) stated by (VII-8). 

3 

([.] means largest integer, [5 •£■] = 5). 

The operators W and W are identical to the operators T(x) and T(x*) as 

A 

defined by (VII-7) respectively. W and W are closed operators. The properties 
of the resolvent of W are indicated in (iii) and (iv). 

In the case where x = x* a more explicit result follows from the 

theorem: 

Theorem VII-4. Let t be a strongly elliptic formal partial differential opera- 
tor as given by (VII-1) and satisfying (VII-2) . Let ft be a bounded subdomain 

ft C ftp. Let 3ft be a smooth surface and let no point in 3ft be interior to the 

2 

closure of ft and let x = x*. Let T be the operator in Hilbert space L (ft) 
defined by 

P(T) = {f e C°°(ft) | f(x) = 3^(3 ft) f(x) = ... = 

( 8 a) f(x) =0, x e (VII-9) 

Tf = xf , f eV(T). 
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Let W be the closure of T, then 

i. The operator W is self-adjoint, 

ii. The spectrum a(W) is a sequence of points {X^} tending to «>, and 
for X in p(W), R(X; W) is a compact operator, 

iii. The operator W has a complete countable set {<t> n } of eigenfunctions. 
Each eigenfunction satisfies the partial differential equations 
t ij> n = ^ > has infinitely many continuous derivatives in 

the closure of the domain ft, and satisfies the boundary conditions 
defining t?(T) of (VII-9). 

The results of this theorem are important for the wave equation problem 
of Chapter VIII. 

B. Stability of the Solutions to an Evolution Equation with 
Strongly Elliptic Partial Differential Operator 

The Dirichlet boundary conditions appear in a large class of boundary 
value problems for elliptic partial differential operators. In stability 
investigations of systems represented by partial differential equations, the 
stability analysis is generally carried out for systems perturbed from 
equilibrium. In this case one can often introduce zero boundary conditions. 

The formulation of the general conditions for the stability and 
existence of solutions to the class of strongly elliptic partial differential 
operators will cover a large class of systems. 

Consider the system as formulated by (II-l) , thus 

~=Lu=-xu (VII-10) 

a t 

where t is now a time- invariant strongly elliptic formal partial differential 
operator of even order 2p with real coefficients defined in a domain ty) of R n , 
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t - J a (x) 3 J . (VII-11) 

|J|<2 P J_ 

Then in any bounded open domain ft, such that ft C ft^, Garding's Inequality holds 
for all u e C^Cft): 

< tu,u>q + K< u,u>q > k < u »u>p (VII-12) 

for some K < and k > 0. 

Next we want to associate with L and thus t a Dirichlet boundary value 
problem. Hence let ft be a bounded subdomain of ftp with ft interior to ftp. Let 
3 ft be a smooth surface with no point of 3 ft interior to ft. Then (VII-10) 
defines an operator differential equation 

= - T u , u e P(T) (VII-13) 

at 

0 (T) = {u e c"(ft) | u = 3 ( 3 ft)u = ... = ft) u » 0; x e 3 ft} (VII-14) 

T u = tu ; u e P(T). 

The stability theory developed in Chapter V cannot be applied to 
(VII-13) even though it is an operator differential equation. In order to 
get (V-l) T must be defined on a complete space with suitable norm. Introduc- 
tion of distributions and Theorems VII-3 provides the appropriate differential 

2 

operator W, the closure of T with respect to L (ft) , as 

W f = Tf ; f e P(W) 

P(W) = Hg(ft) 0 H 2P (ft) C L 2 (ft) . (VII- 15) 

And the corresponding operator differential equation becomes 

— = - W f f e P(W). (VII-16) 
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Furthermore R(W) C L^(fi) and since is dense in L^(ft) Z?(W) = L^(J?)» thus 

(VII-16) is an operator differential equation, associated with (VII-10), which 
belongs to the class (V-l) . 

Theorem VII-3 enables us to define the adjoint of W, W*, as the 

A A _ ^ A A A 

closure W of T where i/(T) = V(T) and T u = x* u for u e P(T). Since x is 

strongly elliptic x* will also be strongly elliptic, thus satisfying Garding's 

00 

Inequality for u e C Q (Q). Replacing x by x* in Theorem VII-3 gives 
t?(W) = P(W*) = H^(Q) = t?(W) . This is an important property, to be used 

in the following stability theorem. Notice that if P(W) = V(\ -7*), then W is 
self-adjoint if x = t* 

It should also be noted that the imbedding of the closed subspace 
in L (fi) implies that there exists a constant C > 1 such that 

< f ,f> p > C < f,f> Q for all f e H^Cn). (VII-17) 

Theorem VII-5 . Let x be a strongly elliptic formal partial differential oper- 
ator as defined by (VII-11), satisfying the Garding Inequality (VII-12) . Let 
T of (VII-14) be the operator associating x with a Dirichlet boundary value 
problem, and let W be the closure of T as defined by (VII- 15), 

Then the null solution of 

= - W f f £ V (W) — H 2p (fi) (VII-18) 

dt u 

2 

is asymptotically stable with respect to the L -norm if a > 1 can be found 
such that 

i. I |f I I* > C 0 ! Jf | !g feV(W) (VI 1-19) 

and 

ii. C Q k - K > 0 (VII-20) 
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where k and K are two constants satisfying Garding's Inequality for x. 

Proof : Define the symmetric, bilinear, positive definite functional V(f,g) by 


V(f,g) = < f»g> 0 for f »8 e h 2 (fi). 

The Lyapunov Functional v(f) is thus defined by 

v(f) = V (f ,f ) = < f,f> Q ¥ f e L 2 (fi). 

Then v(f) = 2V(-Wf,f) = 2< -Wf,f> Q f e P(W). (VII-21) 

Since x satisfies (VII-12) we have by Corollary VII-1: 

< Wf,f> Q + K< £*£>0 > k < f,f> p V f e V(W) 

and with (VII-19) : 

< Wf,f> Q > (C Q k - K) < f,f> Q V f e V(W) (VII-22) 

W is closed, and so is -W, P(W) = P(-W) = L 2 (n). 

Since P(W) = P(W*), it follows from (VII-22): 

< W*f,f> Q = < f,Wf> Q = < Wf ,f> Q > (C Q k - K) < f,f> Q (VII-23) 

for all f e P(W*) = V(W), For (C^k - K) > 0»W is a one-to-one mapping from 
Hq(&) f) H 2p (Q) into L 2 (ft) and R(W) is closed. Hq(£ 2 )H H 2 ^(ft) is reflexive ^ 2 ^ , 

thus W* maps H 2 ^(S2) into L 2 (ft) and since N(W*) = R(W) J ’ , the annihilator 

9 J. 

of R(W) in L^ft), it follows from the closedness of R(W) that R(W) = N(W*) . 

However, by (VII-23) also W* is one-to-one and W(W*) - {0} , thus R(W) = L 2 (ft). 
Substitution of (VII-23) into (VII-21) gives 

❖(f) = 2 < -Wf,f> Q < -2(C Q k - K) < f,f> Q Vf e V(W). 

(VII-24) 

— — o 

Thus for Cpk-K > 0 the operator -W with 'O(-W) = L (ft) is strictly dissipative 
2 

and R(I+W) = L (J2). Hence by Theorem V-5 the Lyapunov Functional v(f)=< f,f>Q 
assures the asymptotic stability of the null solution of (VII-18) . Moreover, 
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-W generates a semi-group {T^ 


t > 0} in L^(fl) with 


I Mo =< e ' at t>o, 

where a = C^k-K. 

In the stability analysis the objective becomes thus to reduce <-Wf,f>Q 

to the form - o< f,f>Q. The choice of the maximal Cq which satisfies (VII-19) 

(31) 

is very important. Integral inequalities are available to facilitate this 

reduction. These integral inequalities are to some extent based on an estimate 

of the eigenvalues. In Chapter IX it will be shown that a proper selection of 

v(f) can improve considerably the effectiveness of the use of these inequalities. 

A special case arises when x = x*. Then it follows from Theorem VII-4 

that the spectrum of W, uniquely determined since W = W*, is a sequence of 

points (X } tending to °°. Let X . be the smallest X ; . then 
n min n 


X min < f,f> 0 < < Wf,f> 0 V f e P(W). 

Thus, if X > 0 we get 


< - Wf ,f> n < - X . < f ,f> n <0 V f e D(W) 

U = min u 

and -W is strictly dissipative. However, the determination of the eigenvalues 
is generally not easy, especially when the coefficients of x depend on the 
system parameters. In many such cases integral inequalities might be used more 
easily. 

The question can be raised if similar results can be derived for the 

asymptotic stability of the null solution to (VII-16) with respect to the H 01 - 

norm. This question becomes important when the Sobolev Imbedding Theorem is 

used to deduce stability properties with respect to a C^-norm. 

2 

From Theorem VII-3 follows that the closure of T in II is VJ„ . For 
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2 

Garding's Inequality to hold in H (ft) it appears that P(W 2 ) must be 
restricted to 

P(w 2m ) = HP +m (fi)fl H 2p+2m (fi)C L 2 (ft) 
and giving for all f e P(W 2m ) 

< W 2m f - f> 2m + K< £ - £> 2m S k < *’ £> 2m+ p- < VII ' 25 > 

This gives rise to the following conjecture: 

Conjecture VII-1 . Let W 2 be the closure of T, as defined by (VII-14), with 
2 

respect to H (ft) . Then the null solution of 


tI * - W 2m f £s V( V ~ O®* n " 2p+2m <« 

2 

is asymptotically stable with respect to the H -norm if there exists a 
C 2m > 1, satisfying 

£ - H f ll^2 n i C 2mll £ llL V££P < H 2») 

and 


ii. C. k - K > 0 

/m 

where k and K are the constants satisfying Garding's Inequality for t. 

The derivation of this result is analogous to the proof of Theorem 

VII-5. 

Then for sufficiently large 2m > t^] + j + 1 it follows from Sobolev's 
Imbedding Theorem that f e H 2m (ft) implies f cC^ (ft) and 



( 0 ) 


< C | | f | 


2m 


where c is a positive constant. 
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However, for some initial function f^e m ) = 

(23) 

remains in ^ or t - » Jit ' ts t * ie Sem:i - - 8 rou P generated by T ^us 

for sufficiently large 2m >[y] - 2p + j + 1 it follows by Sobolev’s Imbedding 
Theorem, that _T f ^ e implies J^f^ e (fi) and 


I ^t f 0^C J (fl) = c ll l-t f ol l2m+2p 

where is some positive constant* Thus any asymptotically stable semi-group 
trajectory in H 2in+2p (Q) should be asymptotically stable in (H) . 

Next consider as a specific example a diffusion equation 


— = - b u ; t > 0 (VII-26) 

9 t 9 x 

where c is a positive constant and b is a constant, 0 < x < 1. Assume the 

boundary conditions to be u(t,0) = u(t,l) = 0. (VII-26) is thus formulated as 

the Dirichlet boundary value problem for a second order partial differential 

equation defined on [0,1] of r\ The results are compared in Chapter VIII 

9 2 

with those obtained for the case where L = -c — „ + b is taken in the wave 

^ 


equation representation. 


Here x = -c — r + b = x* is a strongly elliptic formally self-adjoint 
9x 

partial differential operator defined on [0,1], 

(VII-26) is associated with the closed self-adjoint operator differen- 
tial equation of class (V-l) : 


d f 
dt 


= - W f 


f e P(W) 


V(W) = hJ([0,1]) Cl H 2 ([0,1]) 

W f = xf f e P(W). 


(VII-27) 
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Take as Lyapunov Functional: 

v(f) = < f,f> Q = | | f | | q f e fl(W) 

evaluating v(f) on the c”([0,l]) functions 

12 1 

v(f) = 2 < Lf,f> = 2 /{c|-4 • f - b f 2 >d x = -2 / { c ( |£> 2 + 

u o 3x 0 dx 

+ b f 2 } d x. 

(31) 

This last expression can be reduced with the integral inequality 

r 1 , 3 f N 2 . 2 f 1 ^ 2 . 

) d x > ir z J f dx 

0 dX "0 


v < -2(c tt 2 + b) / f 2 d x = -2(c ir 2 + b) < f ,f> . 
“0 U 

A sufficient condition for the asymptotic stability of the null solu- 


tion of (VIII-27) is thus 


c ir + b > 0, 


Or since for the strong ellipticity of t,c > 0 is required: 

2 

c > 0 and b > - c ir . 

2 

Thus — W generates a semi-group { ; t > 0} in L ([0,1]) which is bounded by: 


lit 1 1 0 ; 


t > 0, 


where 


a = b + it c. 


The results obtained are identical to those obtained by evaluating the 
eigenvalues of -t. These are 


X n = b + cn 2 ir 2 (n = + 1, + 2, ....) 

2 9 f 2 2 

Negative b are allowed as long as b > -c ir , since c ( - — ) + bf =0, implies 

O X 

for f(0) = f(l) = 0 that f = 0 is the only solution. 
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VIII. STABILITY OF WAVE EQUATIONS 


In this chapter the stability properties of an important class of 
partial differential equations will be investigated. Wave equations appear 
frequently in the mathematical representations of physical systems • The 
development of a Lyapunov Functional for this system is sketched in Section A. 
Specific applications are given in Section B. The formulation of the Lyapunov 
Functional and its derivative show that the system possesses the group property 
rather than the semi-group property, and thus necessary and sufficient condi- 
tions for stability can be established. 

A. Development of a Lyapunov Functional 

The follcwing formal derivation of a suitable Lyapunov Functional, 

not mathematically rigorous, constitutes an important step in solving the 

problem of asymptotic stability for the wave equation with Lyapunov’s Direct 

Method. Let us consider the simple wave equation 

2 2 

LJL. + a L“ + b u - c = 0 (VIII-1) 

at: at ax 

with a, b, and c positive constants. Let x be normalized to the interval 
[0,1] and assume u(t,0) = u(t,l) = 0. Under these assumptions one expects 
the null solution u = 0 to be asymptotically stable. However, if the total 
energy of the system is chosen as Lyapunov functional 

v (u) = / {( |-^) 2 + b u 2 + c( ~~) 2 }d x (VIII-2) 

0 at ax 

then, with the assumption that the derivative can formally be determined: 

v(u) = -2a / ( d x. 

o at 

Since v(u) is only negative semi-definite, the null solution u = 0 is assured 
to be stable but not asymptotically stable. 
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A similar result follows if one interprets (VIII-1) in the context of 
(24) 

a dissipative system . Here (VIII-1) is written as 


§T = 


fb u 

f U y 


0 0 >/b 

00 f h 
h -• 


0 0 


Thus the choice of (VIII-2) as Lyapunov Functional is apparently not 

the correct one for showing asymptotic stability. The problem is to select a 

Lyapunov Functional v(u) with a negative definite time derivative. 

The procedure for the construction of such a Lyapunov Functional 

follows to some extent the construction of Lyapunov functions for systems of 

( 2 ) 

linear ordinary differential equations . 

Notice that (VIII-1) can formally be solved by a separation of vari- 

( 22 ) 

ables argument v . Putting u(t,x) = T(t) X(x) (VIII-1) can be replaced by 


the system of equations: 


+ a — + X T = 0 
dt dt 


d 2 y 

c + bX-XX=0 
dx 


and boundary conditions X(0) = X(l) = 0. 

Suppose that the second equation of (VIII-3) is solved and a *| = X^ 
is determined. For this particular example X^ will be real, because the 
differential operator is self-adjoint. Substitution of X^ in the first equa- 
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tion of (VIII-3) gives: 


d 2 T d T , . 

1 r. + a ■ + X. T = 0. 

dt^ dt 

This is a second order linear ordinary differential equation with real 
coefficients that can be written as: 


d T 


1 


dt 



or 


d T, 

£ 

dt 


" - A 1 T 1 - 3 T 2 


d T 
~~dt 


F T 


The construction of a Lyapunov function for this system follows the usual 
procedure. Let 

V(T) = T T P T 
T 

where JP is symmetric and JC denotes the transpose of T then 


d T T T 

— dt — = 1 (£' L + L VI = - t q t . 

The system is asymptotically stable if for a positive definite symmetric mat- 
rix Q, P as a solution to the matrix equation 


F T P + P F = - £ 

is uniquely determined by Q and is positive definite symmetric when is. 

0} is usually taken as the identity matrix. Here the objective is, however, 

to identify X.. with L = - c — _ + b, and v(u) and 'fr(u) should become 

3 x 

quadratic forms which are equivalent to the same norm, or can be reduced to 
equivalent norms. Thus let 


Q 



2 aX 
0 


0 


N 


2aj 
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then P follows as 


2X n + a 4 . 


The next step is to interpret 1? and Q in terms of L. The first 
consequense is to write (VIII-1) as 


d a 
~dt 


- L n 


where 


u = 


> 


/ \ 


"| 

U 1 


u 


0 , 1 




11 

-J 

• 9 \ 


U 2 


u t 


— Ij , —a 

k J 




s J 


Next v(u) must be patterned after V(£) = _T ]? T^ by replacing \ ^ with L and at 

2 

■inner proaucc. xnus 

N 


the same time taking the L -inner product. Thus 


v(u) = / U T 

0 


2L + a 
a 


u dx (VIII-4) 


Integration by parts and substitution of the boundary conditions gives 

^ 9 9 0 9 

v(u) = J {2 c ( r + 2 b u£ + a* U]L + 2 a u^ + 2 upd x. (VIII-5) 


With a, b and c > 0, (VIII-5) is equivalent to the norm of the product space 

d V (T) 


hJ([ 0,1]) x H° ([0,1]), 


For the case of self-adjoint L, v(ju) can be related to 
In this particular example 


<fCu) = -2a / u* 
0 


L 0 
0 1 
S s 


1 3u 


u d x = - 2a / {c (~ -) Z + buif + u^}d 
“ i 3 x 12 


(VIII-6) 


Thus v(u) is equivalent to the norm of the product space Hq([0,1]) x Hq([0,1]), 


83 



For nonself-ad joint operators L (VIII-5) can still be taken as Lyapunov 

functional, but v(u) must formally be evaluated. An example is the panel 

f 16 ^ 

flutter problem analyzed by Parks'" . The Lyapunov functional derived by 
solving a variational problem is exactly equal to (VIII-4) with L replaced by 
the appropriate differential operator* 

The Lyapunov Functional developed above enables the formulation of the 
stability problem for the wave equation in terms of the hypotheses of Theorem 
V-6. Hence we are able to show not only that the solutions are asymptotically 
stable but exist and satisfy the group property. 

Some applications are given in the following section* 

B. Application to a Class of Elliptic Partial 
Differential Operators 

Consider the general system equation: 

2 

— j + a — + L u = 0 (VIII-7) 

at at 

with a constant a > 0 and L = x is a strongly elliptic partial differential 
operator of even order 2p with real coefficients defined and uniformly bounded 
in a domain ft^ of R n : 

T = y a (x) 3 J . (VIII-8) 

Ul,<2p J ~ 

and let x = x*. Then in any open bounded domain 9 such that 9 C 9^ Garding f s 

oo 

Inequality will hold for all u e (ft) : 

< tu,u>q + K<u,u>q > k<u,u>p (VIII-9) 

for some constants K < ® and k > 0. 
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With L is again associated a boundary value problem. Thus ft is a 
bounded subdomain of ftp with ft interior to ft^. Let 3 ft be a smooth surface 
with no point of 3 ft interior to ft. Then x defines a differential operator T 
such that 

P(T) = { u e C°(ft) | u = 3 v .(3 ft) u - ... aj _1 ( 3 ft)u = 0;u e3 ft) 

T u - x u ; u eP(T) (VIII-10) 

T is again not defined on a complete space. In order to get the proper 

differential operator we introduce distributions and let, by Theorem VII-4, W 

2 

be the closure of T with respect to L (ft) so that W is self-adjoint 


P(W) = H P (ft) fi H 2p (ft) 

W f - xf : f e fl(W). 

And the corresponding operator differential equation 


d 2 f . d f . „ . n 

— -r + a + W f = 0 

dt dt 


f e P(W). 


(VIII-11) 


(VIII-12) 


From Corollary VII-1 it follows that for all f eP(W) 


< W f ,f> 0 + K < f ,f> 0 > k < f ,f> . 
Then (VIII-12) can also be written as 

d f 


(VIII-13) 


3 E - w f, vm = H p (ft) fl H 2p (ft) x H p (ft) (VIII- 14) 


where 


f = 


✓ > 


> 


s 

\ 

f l 


f 

; W = 

0 

1 

f 2 

< > 


1-h 

rt 

\ 


-w 

k 

-a 

/ 


(VIII- 15) 
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Since 

then R(W) = It P (ft) x R(W). 

Theorem VIII-1 , Let W(VIII-ll) be the closed self-adjoint differential 
operator associated with the Dirichlet boundary value problem for the strongly 
elliptic formlly self-adjoint x as given by (VIII-8) , then the null solution, 
f^ = _0 , of the system 

£ = W f_ ; V(W) = H P (fi)fl H 2p (fi) x H P (n) (VIII- 16) 

is asymptotically stable if a > 0 and if there exists a > 0 satisfying 

i. I l f il Ip > c 0 ll f illo V P(W) (VIII-17) 

and 


d _f 0 1 

“dt = -W -a 



ii. k - ~ > e > 0 (VIII-18) 

°0 

where e > 0 and k and K are the constants satisfying Garding's Inequality for 
x (VIII-9). W generates a group of exponential type in H P (fi) x H^(n). 


Proof: As the derivation in the preceding section shows, define a bilinear 


functional 

A 


2 2 

»«.£> - l l < q. p ±1 s.i > o 

i=l j=l 1 J u 


f i» e c 0 


with 


P = 


2L + a 


a 

2 


Since the coefficients of L are uniformaly bounded on ft there exists a constant 
k^ > 0 such that by Schwarz’s Inequality: 

I < ^ S^qI = ^1 11^1 lip I I®1 1 Ip ^1* ®l G C 0 
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and 


lW>l;2 1‘ 1 l|i 1 ll p lls 1 ll I , + a 2 ||£ 1 || 0 || gl || 0 + . IlfJ^ ||g 2 || 0 + 

+ a||£ 2 || 0 || gl II 0 «|| f 2 II 0 ||, 2 || 0 

< (2k 2 + a 2 ) | |f x | | p | |g x | | p + a| |f x | | p | |g 2 | | 0 + a| 1 f 2 | | 0 | |g x | | p+ 

+ ? 1 1**1 lo Iloilo 

<2k x ||* x Mp ll*illp + I Mo ll* 2 Mo + fe H f illp + 

+ I l f 2 l l 0 } { a l l g il Ip + I ls 2 l lo } 

< = { ^1 |f x | | p + I |f 2 | | 0 }{ J2k ± I | gl | | p + I |g 2 | | 0 > + 

+ { a ||f 1 M p + ||f 2 |l 0 } {a|| 8l || p + ||g 2 || 0 } 

< [{max( s/Ik^l)} 2 + (max(a,l) } 2 ] [ | | f x | | p + I l f 2 l IqHI Is^ l p + 

I I g 2 1 1 0^ 

where we have used the fact that 

|| fl || p > | | fl | | Q and| | gl | | p > || 8l || 0 . 

Next let k 2 = {max( \/2k^,l) } 2 + {max(a,l) } 2 , and take the square on both 
sides of the inequality sign: 

|v<£,jl)| 2 t||t 1 ll p + IUjIIo )H||g x || p + 1 1 S 2 1 1 0 > 2 

-k 2 (||f 1 || 2 + 2||£ 1 || p ||£ 2 || 0+ ||f 2 || 2 H| |g x | | p + 2 1 |g x | | p ||g 2 || 0 + 

l|g 2 llo 1 
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with the inequality: 


2 1 a. b | < (v a 2 + v - '*' b 2 ) 
which is valid when v > 0, there follows 


|V(f,£)| 2 < k 2 {(v + Dllfjlp + (l+v" 1 )||f 2 Mo H(v+ 1 ) I | g;L | Ip + 


Thus 


+ (i+v h \\g 2 \\ z 0 ) 

< k 2 (max(v+l, 1+v -1 )} 2 | \f_\ | p>() | |g|| p}0 • 
for D > max(v+l, 1+v ; > 0 there follows 


|V«,4)| ,< D||£|l p>0 IWIp.o 

where | | • | | n denotes the norm induced by the inner product: 

D * U 


<£.£^,0 " [ l | J |, p aJ f l ‘ 3J S 1 + f 2 ‘ 8 2 )d x 


for all £ ,j» in C~ (fi) x C~ (fi) . 

a 2 

Furthermore V(£,£) = 2 < f^,Lf^>Q + a < + 

+ 2a< f 1 ,f 2 > 0 + 2< f 2 ,f 2 > 0 1 e c“(0) x C~(fi) . (VIII-19) 

With Garding's Inequality (VII-9): 

V(f,f) > 2k < f 1 ,f 1 > p - 2K< + a 2 < f 1 ,f 1 > Q 


+ 2a < f 1 »f 2 > 0 + 2< f 2 ,f 2 > 0 — e C 0 ™ x 

And with (VIII-17) and (VIII-18) 

v (£»£) > 2 e < f 1 ,f 1 > p + a 2 < f 1 ,f 1 > 0 + 2a < f 1 »f 2 > 0 

+ 2< f 2 ,f 2 > o L e x <£(n>. 
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Thus with e > 0 and a > 0 there exists some constant d > 0 such that 

V(f,f) > d| |f | | 2 Q for all f e C~(fl) x c“(fl). 

By continuity V(f,j») can be extended to all of H P (S2) x Hg(ft) = H(Q) . 

|V(f,^)| = |V(f,£)| <n||f|| p(0 iyi p>0 f,fie H(fl) (VIII-20) 

and V(f,f) > d | |f | | 2 f e H(fl) (VIII-21) 

— p,u 

Then by the Lax-Milgram Theorem there exists a symmetric, positive definite, 
bounded operator _S eL(H,H) such that 

V(f,_g) = s <S f, £> r f., £ e H(fi). 

For functions .f ,£ e V[W) C H, since W is closed extension of L in H 
with domain P(W) = H p (fl)/1 H 2p (fl) it follows from (VIII- 19) that 

v (f»£) “ < f !» (2W + a 2 ) g x >o + < f lt a g 2 > Q + 


+ < f 2» a s l > 0 + < f 2* 2g 2 > 0* 
Let v(f) = V(f,f) = < jf ,.S f >h for f. e H(fi) 


(VIII-22) 

(VIII-23) 


Next we evaluate v(f) = 2V(W_f, < f) for f_ e P(W). Notice that actually 
(VIII-23) must be taken for V(f.,f). v(£) can be taken as (VIII-22), if the 

proper restrictions are placed on W, so that V(W f_,f) = V(£, W _f). Taking 
first 


V(W f,f) = 2 < f 2 , W f 1 > Q + aT< f 2t f 1 > 0 + a < f 2 ,f 2 > 0 - 
- a < W f 1 ,f 1 > 0 - a < ^ 2 * f l > 0 2< W f l» f 2 > 0 “ 2a< f 2* f 2 > 0 = 


= -a < W fpf-^g - a < f 0 ,f 9 > 


2*2 0 


f E P(W), 


Next take 
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V(f, W f) = -a< f lt W f 1>0 + 2 < f r W f 2 > 0 + a 2 < f^f^ - 

a < f l» f 2 > 0 2 < f 2» W f l > 0 + a< f 2 ,f 2 > 0 2a < f 2’ f 2 > 0 = 

= -a< f 1 ,K f 1 > 0 - a < f 2 » f 2 > 0 f, e V[U) 

if W is self-adjoint. Since 

V Wp : * w *i. £ i"o s E < f i-V P h z v < w > 

then, for any a > 0, 

-e | |f | |^ „ < v(f) = 2V (10 f ,f) < — E | |f | I 2 f eP(W) 

p,u ” P» U 

for some constants e and E, e > E > 0. Hence there exist some constants 
e E 

a = -j > 0 and 3=^->0, a>8>0 such that 

- aV(f:,f) < v(f) < - 6 V(f,f) f_ e 1?(WJ (VIII-24) 

In order to apply Theorem V-6, there remains to show that R(A _I - W) = Hq(JJ) 
x Hp(ft) for real X and |x| sufficiently large. Take any vector 

£ = col(g^,g 2 ) in x Hq(J 2), then we must show that for |x| large there 

exists an f e Hg(fi) fl H 2 ^(fi) x Hg(ft) such that (X _I - (£) f_ = j». Thus let: 

X f i - f 2 = g i £ Hg(fl) C H°(fi) 

W f ] _ + (X + a) f 2 = g 2 e H°(n) 

since f 2 e Hq(Q), f 2 e Hq(J 2), thus substitution of f 2 in the second equation 
gives 

w^ + x (X + a) f^ = (X + a) g^ + g 2 = g e Hq(J2). 

For | X | sufficiently large X(X + a) is always positive for X > 0 and X <0. 

In Theorem VII-5, it has been proven that R(W) = Hq(&). From the statements 
of the theorem it follows that R(W + X(X + a) I) = for |x| sufficiently 

large. And thus R(X JC - W) = H^(n) x Hq(J 2) for |x| sufficiently large. Thus 
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by Theorem V-6 W generates a group of exponential type (T^j t e (-»,« )} in 
H^(0) x Hq(J 2) and there exist constants °°>M>l>m>0 such that 


« lp,0 5 I lit il lp,0 i M ''"I HI lp,o 


-et. 


And the conditions for asymptotic stability of the null solution f_ = 0_ of 

(VIII- 16) are thus a > 0 and k - yr > e > 0 for some small e > 0 where k and 

L 0 

K are the constants from Garding’s Inequality for t, and Cq is a constant 
given by the estimate 


H f iMh c oll f illo 


£ x e fl(W). 


This completes the proof of the theorem. Thus the main computational difficul- 
ty is find maximum values of k, Cq and minimum values of K. 

Next consider the wave equation 


2 2 

3 u 3 u 3 u , _ 

— r + a — - c — r + b u = 0 

3 t 3 t 3 "xT 


(VIII-25) 


with a, b and c constants, 0 _< x < 1. And assume boundary conditions 
u(t,0) = u(t,l) = 0. The corresponding representation of (VIII-25) in terms 
of a closed operator W gives 

d f 




with 



f f > 


✓ \ 


✓ \ 

£ = 

f l 

= 

f 

; w - 

0 1 


f„ 


f 


-W -a 


2 , 


L tj 


/ 


V[w) = Hj([0,l])n H 2 ([0,1]) x H° ([0,1]) and 


3 2 f 


W 


f= ” c_ ~ 2 ~+ b f for f x e t?(W) = hJ([ 0,1]) fi H 2 ([0,1]) . 
3 x 
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By Theorem VIII-1, the conditions for the asymptotic stability of the null 
solution - £ are a > 0 and < W f^,f^>Q = e < f-^’^Vl * Thus evaluating 

1 a 2 2 

< W * l<tl , 0 - / (- C _il f^ + b f^} d x on the C^([0,1]) functions gives 

0 8 x 

1 * 2f l 9 9 

< L f l> Vo = / { C ( "3T" ) + b f l } d x = 

2 

> s < + “ e)if + b - e} < 


where the inequality 


1 3f 

/ ( 8 -x 
0 d 


^ d x > IT 2 



X 


has been used. Thus the null solution to (VIII-25) is asymptotically stable 

for a > 0 and lim {(c ir^ + b) - e(ir^ + 1)} > 0. 

£->0 

2 

Hence for c > 0 and b > - ir c, will generate a group {T ; t e( -°°, »)} of 
the exponential type in Hq(R) x 

The conditions for asymptotic stability are thus identical to those 
found for the diffusion equation in Chapter VII. The solution to the wave 
equation possesses the group property, while the solution to the diffusion 
equation possesses only the semi-group property. 

The eigenvalues X for (VIII-3) were determined in Chapter VII as 

A n = b + cn^ ir^ (n = + 1, + 2, ...). 

Substitution in the first equation of (VIII-3) gives 

+ a ~ T + (b + cn^ ir^) T=0 n = + 1, + 2, ... 
dt dt 
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The eigenvalues y^(n) and f° r t ^ le characteristic equation are thus 

Since the condition for asymptotic stability is Re p(n) < 0, it follows that 

2 

a must satisfy a > 0. For n > 1, b + cti tt 2 > 0 implies for all n, c > 0 and 

b + c it 2 > 0. 

Hence the result: a>0, c>0,b>-cir 2 . This is the same as that just 

found by using Lyapunov stability theory. 
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IX. APPLICATIONS 


In this chapter a number of applications are given. The emphasis in 
these applications is on the formulation of the problem in such a way that the 
formal mathematical operations yield rigorous mathematical results. The first 
example shows that the Lyapunov stability theory of Chapter V is applicable to 
systems of ordinary differential equations. In the next examples the emphasis 
is placed on the selection of the Lyapunov Functional, so that the application 
of well-known integral inequalities gives the maximal parameter range for 
stability of the null solution. Some results from Eckhaus^ are demonstrated 
using the Lyapunov stability theory. The corresponding nonlinear cases are 
given in Chapter X. 

Example IX-1. Consider the system 


*1 X 2 

5 x = A x , 

*2= “ 2x i ** 3x 2 

x e R 2 , A e L(R 2 , R 2 ) 


/ \ 



X, 


0 

1 

1 

X 2 

; a = 

-2 

-3 



S 


(IX-1) 


Obviously the solution x « () is asymptotically stable, because the eigen - 
values of the coefficient matrix are -1 and -2. 

Next introduce a norm: 

. , . ,2 T 2 , 2 

I 1*1 l 0 = < >0 = - - = X 1 + x 2 * 

2 2 

Then A, a bounded linear operator with V (A) = R and R(A) = R , is well 

2 

defined on the Hilbert space Hq = (R ; < .,.>^) and (IX-1) so defined belongs 
to the class (V-l) . 
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Then take as positive definite bilinear functional V(:c,£) = < x,^>q 


for x,^ e R . And the Lyapunov Functional v(x) is: 

v(x) = V(x,x) = < x»x>q 


x e R 2 . 


Thus 


1 2 
j v(x) = V(A x,x) = < A x,x> Q = - x.^ - 3 x 2 


which is not even negative semi-definite for all x s R . Hence the inner 

product < ♦»* > q does not provide us with a suitable Lyapunov Functional. 

2 

Next consider V(x,^) = < :c,£>^ = < _x> J3 ^>q, 2£»X e ^ » 

2 1 


with 


S = 


1 1 

v. / 


_S is bounded, positive definite and symmetric, thus < • »»>^ is equivalent to 
<•>•> 0 . v(sc) = < 2 L»x>^» moreover 


Then 


with 


3 - Js 3 + 7 ? 

2 — < 2 W ?0 = < -'->1 = — 2 < 


v(x) = < A x, x>^ - -2x^ - 3 x ^2 - 2x 2 


(IX-2) 


~ ~2 K XjX^q < < A jc , —>1 = “ 
Combining (IX-3) with (IX-2) gives 


2 < £*£>0 ‘ 


(IX-3) 


- 9. .2 V(x,jc) < v(x) = 2V(A x, jc) < - .191 V (x , jx ) . 

2 

Since R(I - A) = R , A generates a group of exponential type 

2 

{T fc ; t e(-~, »)} in = (R , < . ,.>^) such that 




-,0955t 


£l I ]_ 
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And the origin it = £, is asymptotically stable. The selection of the bilinear 
functional V(jc,^) is crucial here. 

The following two examples were studied by Eckhaus ^ using approxi- 
mate methods. The linear cases are given in this chapter, the nonlinear cases 
are analyzed in Chapter X. 

Example IX-2 . Burgers* model to describe turbulence as studied by Eckhaus 



(IX-4) 


0 < x < 1, R constant and boundary conditions u^(t,0) = u^(t,l) = 0. 

The linearized system is: 

9 u 

~ =L £ (IX-5) 

where 

u = 

0 < x < 1 and boundary conditions u^(t,0) = u^(t,l) = 0. The operator 
1 9 ^ 

t = - 1 - — — zr in Lis a formally self-adjoint strongly elliptic partial 
K _ Z 
3 x 

differential operator for R > 0. In order to define ^ on a complete space we 

must introduce distributions. Notice that u^ of (IX-4) is independent of x f 

thus L is defined in the Cq([0,1]) x R functions. to be introduced next, 

2 1 

will thus be defined on the Hilbert space L ([0,1]) x (R ; <*»*>q), where 

2 2 
< u,u>q = u . This space will formally be indicated by L ([0,1]) ; I ) . 
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Following Chapter VII, we get for (IX-5) 


with 


W = 


d JL 

d“t 


W 0 


0 k 


= - W £ 


r 


-L 

0 


1 

R 


for feV (W)= hJ ([0,1] ) f) H 2 ([0,1]) x (R 1 ; <*»*> 0 ). W is 

operator with P(W) £ L 2 ([0,1] ; I) and R(WJ = L 2 ([0,1] j _I ) . 

belongs to the class (V-l). 

Take as positive definite bilinear functional V0£,j»): 

1 

v “ I dx + f 2 8 2 — >£ E 1 

The Lyapunov Functional v(f ) becomes thus 

v(f) - V(f,f) = < f ,f> Q „ f e V[W) 

The next step is to evaluate v(f) = 2 < - (l[ .f ,£ >q °n the Cq ( [0, 

* — ~ f f-2 , 1 p 3 ^f*i i j 1 .2 

< L - , - > 0 ~ l {f l + R f lT T~ X “ R f 2 “ 

0 3 x 


3 f 


_ , r f 1 , °%2 .2, , . 1 .2, 

” [. R 3x - f i> d x + r f 2 ] . 


0 
(31) 

The integral inequality 

1 3f, , 1 9 

/ (“air ) 2 d x > ir 2 / f d X 
0 3X 0 1 

holds for the first term in the integrand, thus 

1 .2 


L f,f > 0 <-[/{—-!) f \H x + | f 2 ] 


0 

Thus - W is strictly dissipative for 


0 < R < TT^ 


(IX-6) 


a closed 
Hence (IX-6) 

. 2 ([0,1] ; I). 

,1]) functions: 


(IX-7) 
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and R(JE + W) =L 2 ([0,1] ; I_), hence by Theorem V-5 the null solution to (IX-6) 
is asymptotically stable. W generates a contraction semi-group (T t ;t> 0} in 
L 2 ([0,1] ; I) with 


IMo* 


-at 


t > 0, 


v 2 1 

a = min ( •= 1, ■=■) . This condition is identical to that obtained by 

R R 


requiring that of L be positive, since from Eckhaus 

r2 


Hn " S’ (n+1) - 1 


( 6 ) 

n ” 0 ) 1) 2 j ••• 


X 20 R s 


Example IX-3 « Next consider a second example from Eckhaus 


( 6 ). 


3 U 


3 u * 2 2 v 2 v ^ -l. 

9t" (x + ^ )U " 7^ X 3“"R 


1 3 2 u , 1 3 u 


3 x 


+ — - + 

2 2 3 x 


+ R 2 [ j u 2 d x] u = 0 (IX-8) 

0 

0 < x < 1, R a positive constant, and boundary conditions u(t,0) = u(t,l) = 0. 
The linearized system is given by 


3 u 
3 t 


L u 


(x + 


h )u 


+ 


2 

n 


X 


3 u + _1 
3 x R 3 x 


u 


(IX-9) 


0 < x <_ 1 and boundary conditions u(t,0) = u(t,l) =0. x = -L is for R > 0 
a strongly elliptic formal partial differential operator, which is, however, 
not formally self-adjoint. (IX-9) can again be formulated in terms of the 
operator differential equation (V-l) by letting: 

= - W f (IX-10) 

with W f - -L f for f e V(W) = Hj([0,l])n H 2 ([0,1]). Take as positive 

definite bilinear functional V(f,g): 
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V(f,g) = < f,g> n . 


f, g e l/([ 0 ,l]). 


Then we have to evaluate again V (— W f,f) on the Cq([ 0,1]) - functions. 
Integration by parts and making use of the integral inequality (IX-7) gives 

V(L f ,f) < - ( ~ - ^ - x 2 ) < f,f> Q . 

Thus a sufficient condition for the asymptotic stability of the null solution, 
f= 0 ,is 

ir 2 1 , 2, n 

— - ae-"” . lx ’ =• 0 

xe[ 0 ,l] 

or 

0 < R < (1 + 2it 2 - Jl + An 2 ) . 

One can improve considerably on this last condition as an evaluation 
of the eigenvalues of X suggest (see Eckhaus v J ) by observing that x is 
equivalent to the strongly elliptic operator as given by: 

x = ~ — 7 —* (p (x) ) + q(x) 

e w (x) 3 x r 3 x M 


p = exp yfK x 

w = R exp ft x 2 

2 2 

q = - x . 

An inner product <• »*> w which is equivalent to <• »*>q can now be chosen for 
the positive definite bilinear functional V(f,g): 


V(f,g) = < f,w(x) g> Q = < f,g>^ 


f,g e IT([0,1]). 


It can easily be checked that W^, corresponding to T e , is self-adjoint with 

respect to <•»•> . Evaluation of V(-W f,f) = <-W f,f> on the C~([0,1]) - 
r w e e w O ’ 

functions gives 
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X, 


V(-T f ,f) = - / {e^- ' 8 f ' 2 

e 0 


Jr x 2 2 , 2 .-2-, 


Cft> - R 6 (x" +^)f"}d 

1 

2 


The integral inequality (IX-7) can now be applied to e 
to f. Substitution gives: 


i€x 2 


f, rather than 


v(-T e f, f ) .<-(4-75 > <f - £> » • 

The sufficient condition for the asymptotic stability of the null solution to 
(IX- 10) becomes then 

0 <R < ir 4 . 

As should be expected, this condition is identical to that found by evaluating 
the eigenvalues of x. The choice of bilinear functional V(f,g) is thus seen 
to be very important for the determination of a maximal parameter range for 
asymptotic stability. The equivalence of o»*>q and <«».> w indicates that 
the stability properties of (IX-10) are the same for both. The selection of 
the bilinear functional is also important, when dealing with nonlinear systems. 
This will be shown in Chapter X. 

The following example shows the effect of the elimination of the 
highest order odd derivative of x on the parameter range for which the system 
is stable, by introducing a x g . 


Example IX-4 . Consider the system 

3 u _ , 1 3ji 1 


3 2 u,5 3^u, 1 3 u . 1 N T 

■\ — - H 1- — u) = L u. 


3 t 


R 2 3 x 4 RjR 3 x 3 4R 3 x 2 Vr 3 x 4 


(IX-11) 


0 < x < 1, R> 0 and boundary conditions 


u(t,0) = u(t,l) = 


au 

3 x 


x=0 


= 0 . 


X— 1 
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t = -L is a strongly elliptic formal partial differential operator. (IX-11) 
can be formaulated as (V-l) by defining 


i 


~ W f (IX-12) 

at 

with W f = -L f for f e V (W) = Hq([0,1])H II 4 ([0,1]). As V(f,g) can be 
taken: 

V(f,g) = < f,g> Q f,g e L 2 ([0,1]). 

Then an evaluation of V(-W f,f) on the G~([0,1]) -functions results in a 
sufficient condition for asymptotic stability of the solution f=0 of (IX-12): 


A o 

0 < R < • 

However, x is equivalent to 




(pM 


3 2 v , 1 3 2 ,3. 8 , 1 

— J> + r — 7 + "r 1 — + “ 

3x 3x VR 3 x 4 


x 

with p(x) = e 1 

_ -s-n/R x 
w(x) = R Z e z 


x e is also strongly elliptic. The bilinear functional V(f ,g) should now be 
taken as: 

V(f,g) = < f,w(x)g > 0 = < f,g> w f »g e L 2 ([0,1]) 

<*»*>q and <• » »> w are again equivalent. Evaluation of V(-W e f,f) = <-W £ f ,f>^ 
shows that the sufficient condition for asymptotic stability of the null 
solution to (IX-12) can be extended to 


0 < R <~ * 2 . , _ -l 

x 3 f 4^ x 

The integral inequality (IX-7) is here applied to e ^ and e f 

3 f 

rather than to - — and f respectively. 

a X 
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{ 


The above examples illustrate clearly the application of the stability 
theory developed here. The importance of the equivalent inner products is 
demonstrated for the determination of maximal parameter ranges for asymptotic 

stability of the null solutions. In each case the stability problem has been 
formulated in such a way that all mathematical operations can be formally 
carried out and these formal operations are rigorously justified. 
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X. SOME NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS 


The important advantage of Lyapunov’s Direct Method over approximate 
methods in the stability analysis of finite dimensional systems of nonlinear 
ordinary differential equations is that rionlinearities can be introduced eas- 
ily and without lack of mathematical rigor. In applications, the main 
emphasis is on the so-called "Lur’e type" nonlinearities, frequently 
encountered in control system applications. Similar nonlinearities are found 
in such distributed parameter systems as nuclear reactor systems and quantum 
physics. 

Here again the main problem is that of the existence and uniqueness 
of the solutions, thus a problem in the theory of partial differential equa- 
tions itself. In the following examples, the existence of solutions is not 
rigorously established; instead the following assumptions are made: 

1. The solutions to the linearized system exist and belong to a 
Hilbert space H. 

2. The linearized system is asymptotically stable. 

3. The solutions to the nonlinear system exist and belong for all 


t > 0 to H. 

The sufficient conditions for stability can then be based on estimates 

for the nonlinear terms in the derivative of the Lyapunov Functional. 

Example X-l . Consider first the non-linear case of Example IX-2 and suppress 

(6) ^ ^2 

(as in Eckhaus v J ) the ? — — term of the second equation of (IX-4). This gives: 

9 , N 


3 u .13 2 u 3u 2 

— = u + - — R 

3 t 3 x 3 x 


/ u d x 
0 


u = L u 


(X-l) 


0 < x < 1 and boundary conditions u(t,0) = u(t,l) = 0. And let R be a posi- 
tive constant. The linearized system is asymptotically stable for 
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0 < R < IT 2 . 

Define (X-l) again on a complete space, which is assumed to be the space for 
the linearized system, then 


d f 
dt 


- W f 


with W f - - L f for f sufficiently smooth. The domain of the linearized 
differential operator is taken as P(W) = Hq([0,1])0 H 2 ([0,1]), 

Take as bilinear functional V(f,g): 

V (f »g) = < f,g> 0 


Then there remains to evaluate v(f) = 2< L f,f>Q on the C n ([0,l]) functions, 


f,g e L ([0,1]) . 

00 

V 

1 


thus 

1 2 2 

< L f ,f> « - / {-f 2 - if 9 —J + f — + R( / f 2 d x)f 2 }d x. f e P(W) . 

U 0 3 x 9 x 0 

Integration by parts and substitution of the boundary conditions gives: 

< L f,f> n = - /{- f 2 + i + R ( f f 2 d x)f 2 } d x f e P(W). 

0 r\ K d X n 


0 


0 


Application of the integral inequality (IX-7) , the fact that J f dx > 0 for 

0 

f e P(W ) and R > 0 gives: 


< L f,f> 0 < “ ( — ~ 1) < f,f> 


R 


f e VOU), 


Thus the modified nonlinear system has, under the stated assumptions, an 
asymptotically stable null solution for 

0 < R < 7T 2 . 

This verifies the results obtained by Eckhaus^^. 

r 33} ^ u 2 

In it is shovm that a similar result is obtained if the 

is not suppressed. 

Example X-2 . Next consider the nonlinear case of Example IX-3. The system 
is given by 
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9 u ,22* ,2 9 u , 1 9 2 u 1 9 u 2 2 , r 1 2 , , _ ,„ 0 * 

— = (x + — ) u + — x — + — — jr - -sr R [I u d x] u = L u (X-2) 

it Vr /r 3x u !< 2 2 ax n 


0 _< x < 1, boundary conditions u(t,0) = u(t,l) = 0 and R a positive constant. 

4 

The linearized system is asymptotically stable for 0 < R < it . 

The next step is to introduce a formulation of (X-2) on a complete 

space: 


d f 
dt 


- - W f 


(X-3) 


with W f = - L f for f sufficiently smooth. The domain of the linearized 
differential operator is taken as P(W) = 1 Iq([0,1])H K 2 ([0,1]). 

The bilinear functional V(f,g) will be taken as 


V(f,g) = < f,w(x) g> Q = < f,g>. 


w 


f,g e u ([0,1]) 


with w(x) = R exp 'T r x 2 . For the linear part of W we must substitute the 
equivalent differential operator 


T - - -4-s I- (p(x) + q(x) 

e w(x) 5x M w 3x ‘ 


where 


p - exp v/R yS 

-2 2 

q = JR - x . 


Evaluation of V(L.f,f) on the Cq([ 0,1]) -functions gives 

\ Vr x 2 9 f 2 x , 2 2 * ^2 

V(L f,f) = - J{ e ( ~-) 2 - R e (X W ‘ + 

o ax 

, 1 _ >/R X 2 . 9 f 2 . „3 n^R X 2 , r 1 ,2 . , -2 . 

+ — R e f + R e [ J f d x] f dx. 

0 


9 x 


Integration by parts, substitution of the boundary conditions and the 

1 2 

inequalitv (IX-7) together with the fact that / f d x > 0 and R > 0 give 

1 o 

V (L f ,f) < - / [ IT “ 7ft “ xf] f 2 w(x) d x. 
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Since 0 < x < 1, L will certainly be strictly dissipative with respect to 


<•»•> for 
w 


h-* R |f| 


0 


R VR 3 

for all x e[0,l]. Thus the null solution of (X-3) will be asymptotically 
stable for all initial distrubances bounded by 


3 1 7T ^ 

max | f | < y -g- ( ^ - 1) for all f sufficiently smooth. 
xe[0,l] 


/ £ \ 

This result is again similar to that obtained by Eckhaus' , however, 
the use of Lyapunov stability theory enables one to obtain it in a straight- 
forward way, without making many complicated calculations as is the case when 
using asymptotic expansions. 

The assumptions made concerning the existence of solutions are not 

more stringent than those made as justification for the use of approximate 

methods. However, the necessity of a self-contained Lyapunov stability theory 

for certain classes of nonlinear partial differential equations remains. To 

( 22 ) 

what extent the theory of monotone operators enables an extension of the 
stability theory, developed here for dissipative operators, to certain classes 
of nonlinear operators is left as a suggestion for further research. 


106 



XI. CONCLUSIONS AND SUGGESTED FURTHER RESEARCH 


A. Conclusions 

The objective of this research is to establish a Lyapunov stability 
theory for solutions to partial differential equations on a mathematically 
rigorous basis. Because of the complex mathematical nature of partial 
differential equations, a type of synthesis method has been developed. 

Lyapunov's Direct Method for ordinary differential equations is 
generalized to a class of operator differential equations. This generaliza- 
tion is based on the fundamental solution structure of groups and semi-groups, 
which is not restricted to ordinary differential equations. 

Once the stability properties are formulated in terms of the group 
and semi-group structure, it becomes possible to impose the conditions for 
their stability on the corresponding infinitesimal generators. This enables 
the formulation of a Lyapunov stability theory, analogous to Lyapunov's 
Direct Method for ordinary differential equations, for a large class of 
operator differential equations — specifically, for those bounded and unbounded 
operators which are the infinitesimal generators of contraction groups and 
semi-groups. The conditions for stability and asymptotic stability for this 
class of operators is based on the Hilbert space theory of dissipative opera- 
tors. 

The Hilbert space structure, associated with dissipative operators led 
to the pivotal notion of equivalent inner products, which enables one to 
relate dissipativity in a satisfactory way to stability, a norm property. 
Sufficient conditions are established for the stability and asymptotic 
stability of semi-groups and for the null solution of the corresponding 
operator differential equations. For groups this is extended to necessary and 
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sufficient conditions for stability and asymptotic stability, respectively. 

The next development is that of defining partial differential equa- 
tions in terms of the general operator differential equation. For this pur- 
pose a formal linear partial differential operator is introduced. The formal 
partial differential operator is extended to an operator defined in a complete 
space by introducing distributions. This is followed by associating the 
partial differential equation x^ith a boundary value problem. In order to 
maintain the general operator differential equation structure this is only 
carried out for the class of strongly elliptic partial differential operators 
satisfying the Dirichlet boundary conditions. Subsequently, a formulation of 
a class of evolution equations and a class of wave equations is obtained in 
terms of the general operator differential equation. 

These two classes of equations are very important in physics and 
engineering applications, and the formulation as operator differential equa- 
tions enables us to apply the developed Lyapunov stability theory for a 
rigorous determination of their stability properties. 

The relation between the stability theory and the group and semi- 
group structure automatically gives the existence and uniqueness properties 
for the solutions of these equations. For the evolution equation a general 
Lyapunov Functional is formulated in terms of the inner product of the Hilbert 
space on which the operator is defined. Subsequent applications clearly 
reveal the importance of the notion of equivalent inner products and its 
consequences for obtaining maximal parameter ranges for stability. 

For a class of wave equations a Lyapunov Functional is developed, 
which not only clearly gives conditions for stability and asymptotic 
stability of the null solution, but also exhibits the group structure of the 
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solutions. 


In a final chapter the stability properties of some nonlinear systems 
are formally investigated. Except for these last results, the emphasis is on 
a mathematically rigorous approach to the important problem of stability of 
solutions to partial differential equations. That such an approach gives 
only results for a limited class of partial differential equations is not 
surprising. The theory of partial differential equations itself is a field 
of extensive research in mathematics and it is forseen that many of its 
developments might find applications in stability studies. Some suggestions 
for further research are given in the next section. 

B. Suggested Further Research 

Future research in the stability properties of solutions to partial 
differential equations can progress in many directions. Along the lines of 
the stability theory developed in this research it must be pointed out that 
only a small class of linear time-invariant partial differential equations 
have been formulated in terms of the general operator differen- 
tial equation. The possibilities of extending this class must be explored, 
even though the ri^id structure might have its limitations. 

There is also a need for investigating the possibilities of extending 
the developed stability theory to those linear time-varying operator differen- 
tial equations, which might possibly generate two parameter groups and semi- 
groups. 

The potential use of Lyapunov stability theory to establish rigorously 
the stability conditions of nonlinear systems justifies a continued effort to 
extend the results for linear systems. The natural extension of linear 
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disspative operators to the nonlinear case seem to be the monotone operators (^2) 
The suggestion is that there might exist a nautral extension of the stability 
theory for linear dissipative operators to one for monotone operators. 

The stability theory developed is based on the Hilbert space structure. 
The stability results are all with respect to norms induced by an inner pro- 
duct. This provides no limitation for finite dimensional systems, since all 
norms are equivalant. The question must then be raised as to how far this 
equivalence of norm principle can be carried through for infinite dimensional 
systems. In other words, which other norms are equivalent to the one induced 
by an inner product? Is there a Banach space structure other than the 
Hilbert spaces that encompasses more equivalent norms and hence a broader 
stability analysis? 

The stability properties obtained are generally those with respect to 

2 

the L -norm. The extent to which this kind of stability implies stability in 
the classical sense, i.e., with respect to a C p -norm should be investigated. 

This is the direction indicated by Conjecture VII-1. Sobolev's Imbedding 
Theorems will play an important role in such an evaluation. 

Even though partial differential equations appear in many engineering 
applications it is plain that the suggested research problems are very 
mathematically orientated. Unfortunately, the complexity of partial differen- 
tial equations seems to require a highly specialized mathematical background 
in functional analysis, topology, generalized functions, etc. It is hoped 
that this will not deter others from investigating the very important problem 
of the stability of solutions to partial differential equations in the 
future. 
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ABSTRACT 


In this report, necessary and sufficient conditions are obtained for 
a closed linear operator A to generate a group { T fc ; - <» < t < »} in a 
Banach space such that { T fc ; t > 0} is a negative contractive semi-group 
with respect to an equivalent norm. These results are refined to the case 
of a group in Hilbert space. 
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'• It is well known that some systems of differential equations, both 
ordinary and partial, can be reduced to the form 

(1) |f “ Ax (xeP(A)) 

where A is a linear operator with domain, t?(A) , and range, R(A), both in a 
real B-space, X. If A is the infinitesimal generator of a semi-group 
{ T fc ; t > 0} C L(X,X) of class (C Q ), a solution to (1) starting at t ® 0 from 
x * XqE??(A) is given by x(t; Xq) = T^Xq for t > 0 with x(0; Xq) = Xq^. If 
A is the infinitesimal generator of a group, then the above solution is valid 
for - » < t < ». 

Criteria in terms of the operator A (in other words, properties of the 
coefficients of the original system of differential equations) which would 
enable one to deduce the existence as well as asymptotic behavior of solutions 
are desirable. One such result already exists for A to be the infinitesimal 
generator of a semi-group (Theorem 1), The object of this report is to refine 
these results to groups of certain types. The development is in terms of 
real B-spaces but can easily be extended to complex spaces. 

DEFINITION 1 

Let {T fc } be a semi-group of bounded linear operators in a B-space. If 
| | T fc | | < N then the semi-group is said to be equibounded ; if N * 1 then {T fc } 
is called a contraction semi-group . If there exist finite M > 0 and g > 0 
such that | |T fc | | < Me” then (T t > is said to be a negative semi-group; if 
M ■ 1 the semi-group is called negative contractive . 

An important property of B-spaces is that a semi-scalar product [.,.] 
can be defined on any B-space such that [x,y] = fy(x) where f^ is a linear 

2 h) 

continuous functional such that fy(y) = | i y 1 | ^ J 
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DEFINITION 2 


Let X be a B-space with norm | | • | | and let [ . , . ] be a semi-scalar product 
on X. Then the semi-scalar product [,,.]^ is said to be equivalent to [.,.] 
on X iff | | • | | ^ aa d 11*11 are equivalent norms on X. 

The class of operators A needed for the statement of our results is 
described in the following definition ^ . 

DEFINITION 3 

Let A be a linear operator with V(A) and R(A) contained in a real B-space X. 
A is called dissipative with respect to the semi-scalar product [ . , . ] if 
[Ax, x] < 0 whenever xeV(A) and strictly dissipative if there exists a y > 0 
such that 

(2) [Ax, x] < - y [x, xl = - y| M 1 2 (xeV(A)) 

The following results are due to Lumer and Phillips 
THEOREM 1 

Let A be a linear operator with V(A) and R(A) contained in a real B-space X 
and such that V(A) is dense in X. Then A generates a contraction semi-group 
in X iff A is dissipative (with respect to any semi-scalar product) and 
R(I - A) = X. 

COROLLARY 

If A is a densely defined closed linear operator such that V(A ) and R(A) are 
contained in a B-space X and if A and its dual operator A* are both dissipa- 
tive, then A generates a contraction semi-group. 

In the course of the proof of Theorem 1 it is deduced that the resol- 
vent operator of A, R(X; A), exists for all X > 0 and satisfies the estimate 
| |R(X; A) | | < From this fact we can deduce a useful criterion for {T^} 
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to be a group and the solution T^Xq of (1) to approach the null solution as 
t -*■ ». We begin by considering semi-groups. 

THEOREM V 

Let A be a linear operator with V(A) and R(A) both in a real B-space X such 
that V(A) is dense in X. Then A generates a negative contractive semi-group 
in X iff A is strictly dissipative and R(I(1 - y) - A) = X where y is the 
constant appearing in (2) . 

PROOF 

Apply Theorem 1 to the dissipative operator B = A + yl, The result follows 
upon noting that B generates {S t ; t £ 0} 'iff S fc = e Yt T t where {T fc } is 
generated by A. 

From the remark preceding Theorem 1’ it follows that R(p; A) exists for 
all p > -y and | |R(p;A) | | < (p + y)**\ The corollary of Theorem 1 can also 
be extended to strictly dissipative operators. 

LEMMA 1 

Let A be the infinitesimal generator of an equibounded (negative) semi-group 
(T t ; t > 0} in a real B-space (X, ||.||^). Then there exists an equivalent 
semi-scalar product, [.,.] inducing an equivalent norm jl*!^ with respect to 
which A is (strictly) dissipative. 

PROOF 

By hypothesis j | T^ | < Me” where a > 0 is a constant and the constant M may 
be taken to be M > 1. If we set S *» e at T then (S ; t > 0 } is a semi-group 

U U u c 

satisfying | | S t | | ^ < M. Let xeX and define | | • | ( ^ by 



In (3 ) f shown that 

(4) 


I |x| l x < ||x|| 2 <M| |x| \ 1 


(xeX) 


(5) l|S t M 2 <1 . 

By the remark following Definition 1, there exists a semi-scalar product 
consistent with |j.j| 2 such that, by Theorems 1 and 1*, 

[ (A + al)x, x] < 0 (xeP(A)). 

This proves Lemma 1 and establishes a starting point for groups in B-spaces, 
and as a special case, Hilbert spaces. 

THEOREM 2 

Let A be a linear operator with T?(A) and R(A) both contained in a real B-space 
( x » 11 * 112 ) such that P(A) is dense in X. Then A generates a group 
{T. ; - °® < t < »} in X such that {T. ; t > 0} is a negative contractive semi- 
group with respect to an equivalent norm | | . | | iff 

(6) - 6 | 1 3t| | 2 = < [Ax, x] < - y| |x| | 2 (xeP(A)) 

where » > 6 > y >0 and is an equivalent semi-scalar product consistent 

with 11*11* and 

(7) R(l(l - y) - A) - X, R(I(1 + <$) + A) - X. 

PROOF 

Suppose that (6) and (7) are valid. Then B = A + yl and C « -A - <Sl are 
dissipative. As in the proof of Theorem 1' it follows that' {f t ; t > 0} is a 
negative contractive semi-group. Moreover R(p; A) exists for all p > -y and 
||R(p; A) I I < (p + y) _1 . Similarly R(pi A) exists for all p <-6 and 
| |R(p; A) I I = <(|p| These last two conditions imply that A is also the 

infinitesimal generator of a group ^ . 
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/*! \ 


Conversely suppose that A generates a group such that { T fc ; t > 0} is 

a negative contractive semi-group with respect to | | *| |> that is, | |t 1 1 

—1 at ^ K 

,(t > 0) where 8 > 0. It is known that for a group | |T~ | | < Me , where M > 1 
and a can be chosen such that a > $ Define S fc = T^e Qt an d define | |.| | 2 

as in (3) but using ||*|| on the right of (3). Then Il s t ll 2 ~ 1 and l!*!^ is 
equivalent to ||*|| and so equivalent to the original norm. It follows from 
(5) that [S fc x - x,x] 2 < 0 and so upon dividing by t and letting t H we have 


[(- A - al)x, x] j < 0* 


This yields an equivalent semi-scalar product and the left side of inequality 
(6) with 6 = a. To show that the right side is also valid consider 

1 2 


( 8 ) 


[T g e 6s x - x,x] < | |T g e Bs x| | 2 | |x| | 2 - .| |x| , 2 . 


Bs 


Next estimate T e x L as follows 

S 1 L 


Since 



! (e -“ ,S ||T s . t e a(s - t) x| | ) . 


= e^"°^ s max(sup ||Te ot x||, ||x|| 
0<t<s 



sup ||T e aT x|| < sup e (a " 3)T ||x|| < e (a " e)s ||x|| 
0<T<S 0<T<S 


it follows that 


I e 8s x||, 
s 


< max 


(I 


(8-a)s 


M l 2 ) < I |x 


Hence (8) yields 


Ac 

' [T g e x— x ,x] 2 < 0 

which in turn implies the right side of (6) with y=8* 
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Finally (7) follows from Theorem 1 applied to the dissipative operators 
T_ t e at (with respect to | | • || 2 > and T fc e 6t (with respect to ||*|| 1 ) 

COROLLARY 

If A is a densely defined closed linear operator such that P(A) and R(A) are 
both contained in a B-space X and if A and its dual operator A’ satisfy (6) , 
then A generates a group such that {T t ; t > 0} is a negative contractive semi- 
group . 

Theorem 2 can be strengthened in Hilbert spaces so that it holds for 

scalar products. The proof is not a straightforward application of Theorem 

2. The difficulty lies in the fact that if (H, (.,.)) is a Hilbert space with 

scalar product (.,.) then H with an equivlanet norm is not necessarily a 

Hilbert space. For example, a Euclidean 2-space (X, ||*||) with 
2 2 2 

| | (x^x^ || = x^ + x^ is a Hilbert space while (X, | | • | | ^) , where 

||(x^, x 2^ I 1 1 = l x il + * s n0t a Gilbert space. 

We require the following lemma. 

LEMMA 2 

Let (H, (.,.)) be a real Hilbert space and (T^; t > 0) a semi-group on H with 
infinitesimal generator A; then 

lim t ^((T x, T x) - (x,x)) ° 2(Ax,x) (xeD(A)). 

t-*0+ C 

PROOF 

The result follows from the identity 

(T fc x, T fc x) - (x,x) = (T fc x, T t x) - (x, T t x) + (x, T fc x) - (x,x) 
and the continuity of the scalar product. 

THEOREM 3 

Theorem 2 is valid for (H, (.,.)) a real Hilbert space and an equivalent 
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2 2 

scalar product with [x,x] = | |x| | and (x,x) = | |x| |^ . ; 

PROOF i 

»■ ■ —»■ ■ l 

i 

It is only necessary to prove that if A generates a group such that (T fc ;t > 0} j 

! 

is a negative semi-group with respect to | |*| |^, then an equation of the form | 

(6) is valid where [.,.] is an equivalent scalar product. Define by | 

oo 

(9) [x, y] ■ /( T.x, Tv) dt 

o Z J 

By hypothesis, | |t | | x < Me“ Pt (t > 0), where g > 0 and M > 1; hence j 

(10) [x,x] = / | |T t x| | ^ dt < / M 2 e” 26t | |x| | J dt = (M 2 /26) | |x| | 2 

o o 

Since {T fc } is a group, there exist constants a > g and 1/k > 1 such that 
| | Tt" 1 1 1 1 < (l/k)e at . ,By using the fact that M^xl^ > ||T t " 1 ||“ 1 | |x| | x it 

I 

follows from (9) that 

(11) [x,x] > / k 2 e -2at | |x| | 2 dt = (k 2 /2a) | |x| | 2 

o 

We leave it to the reader to verify that is a scalar product. The 

equivalence of the two scalar products follows from (10) and (11) . 

To show that an equation of the form (6) is valid we consider 

n n 

[T x, T x] - [x, x] * lim[ / (T T fc x, T g T t x)ds - / (T g x,T s x)ds] 

c n-*° o o 

n+t t 

= lim [ / (T x, T x) ds] - /(T x, T x)ds 

s s s s 

n-^° n o 

t 

a - /(T X, T x)ds, (t > 0). 

s s 
o 

This last equality and Lemma 1 imply that 

(12) 2[Ax, x] « - ||x|| 2 (xeP(A)) j 

I 

Equations (10), (11), and (12) yield (6) with y = g/M 2 and 6 = a/k 2 . i 
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